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PREFACE 

This report represents a part of the Ontario Ministry of the 
Environment's contribution to the International Hydrological Decade 
(IHD) Representative Basin Program. This Program was initiated in 1965 
to study five drainage basins in southern Ontario which were deemed to 
be representative of larger physiographic areas. The "Geology and 
Ground Water Resources of the Wilton Creek IHD Representative Drainage 
Basin" is the second in a series of several such reports. 



G. H. Mills, Director 
Water Resources Branch 
Toronto, January 1977 
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"2,300". 
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ABSTRACT 

This report describes the geology and ground-water resources of 
the Wilton Creek drainage basin. The Wilton Creek drainage basin was 
chosen for study because it is representative of a large area in south- 
eastern Ontario characterized by a limestone plain with sparse to non- 
existent overburden cover. 

The bedrock is comprised of Cambrian and Ordovician sedimentary 
rocks lying non- conformably above Precambrian rocks of Grenville age. 
The sedimentary rocks consist of the conglomerates, sandstones and 
shales of the Basal Group, and the carbonates of the Simcoe Group. The 
overburden is primarily of Pleistocene origin deposited during the last 
Wisconsin ice advance. 

The carbonate bedrock constitutes the major aquifer. The overburden 
is too thin and discontinuous to have anything but a localized affect on 
the hydrogeology . 

The quantity of water available for use is difficult to predict 
because it is dependent on the number of interstitial openings encountered 
during well drilling and their degree of interconnection. Productivities 
of wells finished in different hydrogeologic settings vary significantly. 
In the basin, well productivities tend to be higher for shallow wells 
then deep wells and higher for wells finished in valleys as opposed to 
upland areas. 

For the interval, 1967-68 and 1973-74, the mean precipitation over 
the basin is estimated to be 37.04 inches. The mean total streamflow is 
16.40 inches of which 3.97 inches or 24.2% is estimated to be baseflow. 
Evapotranspiration as calculated by the Munson-Index method is 22.68 
inches and by the indirect or water-balance method, 20.64 inches. 

The chemical quality of the surface water is good. The ground- 
water quality varies with location and depth of sample. Most samples 
are of the calcium-bicarbonate type. With the exceptions of some of the 
deep well water samples, the waters are satisfactory for domestic uses. 
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Figure 1 . Location of the Wilton Crook drainage basin in southern Ontario. 



INTRODUCTION 
SCOPE OF PRESENT INVESTIGATION 

This study of the ground-water potential and the relationships 
amongst the geology, water quality and water availability was undertaken 
in the Wilton Creek drainage basin by the Ontario Ministry of the Environ- 
ment, as part of the International Hydrological Decade (IHD) Representative 
Basin Program. The geologic and ground-water data used in this study 
were compiled from existing reports {both published and unpublished) , 
water-well records, and detailed field mapping undertaken during the 
summer of 1973. The field work included extensive road geology, selected 
traverses, air-photo interpretation and overburden sampling. Stratigraphic 
cross sections were compiled using data from existing sand and gravel 
pits. Hand augering was employed in areas where exposures were lacking. 
In addition, water samples from selected wells in the basin were collected 
in order to determine the chemical composition of the ground water. 

The bedrock geology was compiled using the guidelines suggested by 
B.A. Liberty (pers. comm.). The surficial geology was supplemented with 
information provided by E. Mirynech (unpublished field notes) . 

PREVIOUS GEOLOGICAL INVESTIGATIONS 

The Precambrian stratigraphy of the Canadian Shield to the north 
and underlying the basin is comprised of a variety of lithologic types 
including meta-sedimentary, volcanic and intrusive rocks. The adjacent 
shield area and Frontenac Axis have been mapped in detail by Wynne- 
Edwards (1962, 1967a) and Hewitt (1963). The Precambrian surface is an 
irregularly peneplained terrain with a relief in some areas of 350 feet 
(Liberty, 1971). The irregular relief is transposed under the Paleozoic 
covering, as represented by numerous Precambrian inliers in eastern 

Ontario. 

Overlying the Precambrian rocks non- conformably are the Paleozoic 
rocks of the late Cambrian to Middle Ordovician periods. The Paleozoic 
stratigraphy, as described by Liberty (1967, 1969, 1971), Peterson 
(1969), and Dolar-Mantuani (1975), consists of a sedimentary sequence 
including basal conglomerates, sandstones, shales and carbonaceous 
rocks . 

The overburden of Pleistocene and Recent age is comprised of glacial, 
glaciofluvial and glaciolacustrine sediments deposited during the last 
ice advance (Wisconsin), and alluvium and muck deposits. The area was 
mapped by Coleman (1936), Chapman and Putnam (1966) and on a more detailed 
scale by Mirynech (unpublished field notes) . 
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GEOGRAPHY 
LOCATION 

The Wilton Creek drainage basin lies in southeastern Ontario 
between longitudes 76° 37" and 76° 57' , and latitudes 44° 10' and 44° 
28'. It is situated about 12 miles west of the City of Kingston and is 
readily accessible by highways #2, #38, #133 and #401 and a network of 
county and township roads. Both the Canadian Pacific and Canadian 
National railways traverse the basin. 

The drainage basin is approximately 56.6 sq. miles in area and is 
essentially linear in shape, running northeast to southwest. The maximum 
length of the basin is 23.8 miles and the maximum width is 3.5 miles. 
The position of the basin with respect to southeastern Ontario is shown 
in Figure 1. 

PHYSIOGRAPHY 

The drainage basin lies within a physiographic subdivision referred 
to by Chapman and Putman (1966) as the "Napanee Plain". This area is 
characterized by a flat plain of limestone dipping gently to the south- 
west, resulting in a series of shallow erosional scarps. Relief in the 
area is in part the result of erosion by stream drainage and erosion by 
past glaciations. In the northern portion of the basin, relief is 
approximately 100 feet, while in the south it is about 50 feet. 

Coleman (1936) suggested that jointing controls the orientation of 
topographic features such as escarpments and river valleys in southeastern 
Ontario. This controlling factor is reflected in the general trend of 
Wilton Creek. The fractures in the limestone bedrock have a prevailing 
northeast to southwesterly trend. 

The overburden is generally sparse to non-existent in the basin. 
Relief in the southern portion of the basin and in low-lying areas such 
as the valley proper, has been modified by localized accumulations of 
overburden. 

The major urban centres, in order of size, are: Harrowsmith, Wilton, 
Morven and Violet. A land use study undertaken in the basin in 1969 by 
the Ministry of the Environment, indicates that the primary industry is 
agriculture, with dairy farming comprising the greatest acreage. 

CLIMATE 

The drainage basin lies within Brown, McKay and Chapman's (1968) 
South Slopes climatic regime. The climate is influenced by the moder- 
ating effects of Lake Ontario. As a result, the growing season is 
typically 190-200 days long and extends from the middle of May to late 
November. The mean maximum and mean minimum monthly temperatures and 
mean monthly temperatures for the Hartington area for the interval 1968- 
1973 are presented in Table 1. During this interval, the mean annual 
temperature was 43. 3°F. The lowest mean monthly temperature, 15.0°F 
occurred in January and the maximum mean monthly temperature, 68.4°F, 
occurred in July. The mean monthly precipitation for Hartington and 
Violet from 1968-197 3 is shown in Table 2. The mean annual precipitation 
for that interval was 36.33 and 36.44 inches, respectively. 
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TABLE 1. MEAN MAXIMUM AND MEAN MINIMUM MONTHLY TEMPERATURES (°F) FOR THE HARTINGTON CLIMATOLOGICAL STATION 

FOR THE INTERVAL 1968-1973 

Monthly Temperature Jan. Feb. Mar. Apr. May June July Aug. Sept. Oct. Nov. Dec. Yearly 

Mean Maximum 2 3.6 25.3 35.8 50.2 63.1 72.6 78.6 77.6 70.0 57.1 40.9 27.2 51.8 

Mean Minimum 6.3 8.8 20.5 32.3 43.7 55.1 58.3 57.5 51.6 41.2 29.7 13.0 34.8 

Mean 15.0 17.1 28.1 41.2 53.4 63.8 68.4 67.5 60.8 49.1 35.3 20.1 43.3 



TABLE 2. MEAN MONTHLY PRECIPITATION (INCHES) FOR HARTINGTON AND VIOLET FOR THE INTERVAL 1968-1973 

Location Jan. Feb. Mar. Apr. May June July Aug. Sept. Oct. Nov. Dec. Yearly 

Hartington 1.88 2.86 2.65 2.01 3.60 3.34 3.16 2.75 3.25 2.89 3.71 4.23 36.33 

Violet 2.20 2.71 2.43 2.56 3.30 2.91 2.44 2.63 3.51 3.05 3.65 5.05 36.44 









TABLE 3 COMPARATIVE STRATIGRAPHY FOR THE WILTON CREEK AREA (MODIFIED AFTER LIBERTY. 1969) 
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BEDROCK GEOLOGY 

With the exception of the Potsdam Formation, the entire strati- 
graphic section as summarized by Liberty (1971) is exposed in the Wilton 
Creek drainage basin (Appendix A and Map 1) . The outcrop exposures 
range from Precambrian to Middle Ordovician in age. The Potsdam Form- 
ation is inferred to be present from water-well logs. 

The stratigraphic sequences as interpreted most recently by Kay 
(1960), Liberty (1969, 1971) and Peterson (1969) are presented in 
Table 3. 

PRECAMBRIAN 

The Precambrian bedrock, Grenville in age, has been mapped primarily 
as gneisses, schists, and granites, crystalline limestone and dolomite, 
amphibolite, paragneiss, volcanics and pegmatite. The Precambrian can 
be found exposed only at one location, northeast of Harrowsmith. At 
this locality, an inlier of granite gneiss with an estimated relief of 
75 feet (as determined from adjacent water-well logs) is exposed. The 
Precambrian and Paleozoic contact at this location is shown in Photo la. 
The contact at a location east of Wilton Creek near Inverary Lake is 
shown in Photo lb. 




Photo la. Contact between Shadow Lake Formation and underlying Precambrian 
northeast of Harrowsmith. 

BASAL GROUP 

Lying non- conformably above the Precambrian is a series of beds 
known as the "Basal Group". This Basal Group represents a shallow water 
depositional facies immediately above a period of erosion. The Basal 
Group is divided into two units, the Potsdam and Shadow Lake formations. 
This sedimentary sequence has an overlapping relationship with the 
Precambrian Shield in the area of the Frontenac Axis, as is suggested by 
exposures of Paleozoics on the flanks of Precambrian highs (Figure 2) . 







Photo lb. Angular unconformity between Precambrian and overlying 
Shadow Lake Formation. Photo taken near Inverary Lake 
north of Kingston. 




Photo 2. A road cut exposing the Potsdam Formation east of the 
drainage basin near Brewers Mills. 



Normally present at the contact between the Paleozoic and Precambrian 
rocks is a basal conglomerate composed of rock fragments reflecting the 
underlying bedrock material. This conglomerate is time transgressant , 
occurring as high in the stratigraphic sequence as the Gull River Formation, 
where the older formations are not present. 
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Figure 2. Generalized diagram illustrating the overlapping of the sedimentary sequence on 
Precambrian highs. 

Potsdam Formation 



The Potsdam Formation, whose type-section is in northern New York 
State, is generally considered to be Cambrian in age, although evidence 
for this is lacking in Ontario. The occurrence in Ontario of the fauna 
Lingula acuminata and Ophileta compacta indicates an Ordovician age for 
at least the upper part of the Potsdam (Wilson, 1946 and Liberty, 1971) . 

Liberty (1971) recorded a maximum thickness of 149 feet for the 
Potsdam Formation in the Know 1 ton Lake-Holleford area. This thickness 
is a result of deposition in a Precambrian low, which has been attributed 
to a "meteorite impact crater". The near circular shape of the "crater" 
is evident from air photos of the area. A detailed study of the "crater" 
by St. John (1968) supports the meteorite impact origin. 

The Potsdam Formation outcrops just to the east of the basin, on 
the west shore of Knowlton Lake. It consists of an even- textured, 
medium-grained, white to red quartz sandstone, cemented by hematite and 
calcite. From water-well logs, the average thickness is about 30 feet 
in the basin area. The Potsdam Formation, where it outcrops northeast 
of Kingston, is shown in photos 2 and 3. 




Photo 3. Vertical cylindrical structures in the Potsdam Formation 
at the Park of Pillars. 
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Photo 4. Contact between the Shadow Lake Formation and the overlying 
Gull River Formation near Inverary Lake. The Shadow Lake 
Formation is represented by a black and green shale. Both 
the lower and middle submembers of Member A of the Gull River 
are exposed. 
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Shadow Lake Formation 

The Shadow Lake Formation, first named by Okulitch (1939), was 
redefined by Liberty (1955, 1969) to include only the red, green and 
black shales and arkoses that lie unconformably above the Precambrian 
and the Potsdam Formation. Okulitch had originally included the lower 4- 
5 feet of the Gull River Formation, an argillaceous limestone. 

The Shadow Lake Formation, which has not been differentiated from 
the Potsdam Formation in the water-well logs, can be found draping the 
Precambrian inlier north of Harrowsmith (Photo la) . The sequence exposed 
(about 2 feet) represents a vertical facies change from a basal para- 
conglomerate to an arkose, then to a shale. The conglomerate fragments, 
in addition to an abundance of quartz and jasper pebbles, reflect the 
lithology of the underlying Precambrian surface. The shale, a green and 
red-mottled claystone with a high carbonate content, contains pockets of 
arenaceous material. The uppermost contact is between the black and 
green shale beds of the Shadow Lake Formation containing the ostrocod 
Ctenobolbina clavigera , and the overlying shaly limestone of the Gull 
River Formation, (Photo 4) . This formation appears to be about 20 feet 
thick in the vicinity of the basin. 

SIMCOE GROUP 

The Simcoe Group (Liberty, 1955, 1963, 1969, 1971) is divisible 
into four formations, which in ascending order are: the Gull River, 
Bobcaygeon, Verulam and Lindsay. Only the first three can be found in 
the basin area, the Lindsay outcropping further west. The Simcoe Group 
is approximately 700 feet thick in the Kinston area, while in the Wilton 
Creek area the thickness is only about 400 feet. 

Gull River Formation 

The Gull River Formation was defined by Liberty (1969) to include 
the strata lying between the Shadow Lake Formation and the fine to 
medium-grained limestones of the Bobcaygeon Formation. The Gull River 
was sub-divided by Liberty (1971) into four members (A to D) . In the 
Wilton Creek area however, the repetition of lithology and the diffi- 
culty in correlation requires the grouping of members B, C and D into 
one mappable unit. Peterson (1969) sub-divided the Gull River Formation 
into a "lower limestone" and three separate cycles on the basis of 
lithology and faunal progression (Table 3). Peterson's identification 
guidelines are not incorporated into this study. 

Member A. . . . Member A has been sub-divided into three mappable 
submembers (Liberty, 1971) in the Kingston area. The lower and middle 
submembers are grouped together by the author on the bedrock map (Map 
1) . The lower submember is gray to brown to green, strongly-weathered, 
soft, crystalline to lithographic limestone with thin dolostone beds. 
Liberty (pers. com.) suggested that the texture of these beds, as well 
as the soft beds higher in the sequence, is related to dedolomization. 
The alkali-reactivity of these carbonate rocks and this dedolomization 
phenomenon is described in Dolar-Mantuani (1975) . This submember is 
about 2-5 feet thick in the Kingston area (Photo 4) . 






The middle submember is a brown to dark brown weathering litho- 
graphic limestone, with a distinctive purplish hue on fresh surfaces. 
In some areas, this submember takes on a mottled, digitate appearance. 
Soft dedolmitized beds (0.5-1.5 feet in thickness), which are strongly 
weathered, occur at intervals throughout the section. Mineral molds and 
calcite crystals are abundant, often resulting in a pitted and character- 
istic wormy, weathered surface. In the basin, the maximum thickness of 
the middle submember is approximately 70 feet. The upper contact has 
been drawn at the base of an unfossiliferous black shale phase of the 
upper submember (Liberty - pers. comm. ) , (Photo 5). 




Photo 5. Contact between the middle and upper submember of Member A 
of the Gull River Formation, represented by black shale. 
Photo was taken near Inverary Lake. 

The upper submember is an alternating gray, sublithographic to 
lithographic limestone with brown, soft to hard, fine-grained to finely 
crystalline dolomitic limestone and dolomite beds. Shaly partings and 
horizons are frequent in the sequence, as are calcite vugs, mineral 
molds, limey conglomerates and calcarenites. Dolomitic limestone beds, 
fcund just above the upper contact of this submember, can be used as 
marker beds, where they are traceable over the upper portion of the 
basin. Bedding is relatively uniform and varies from 0.1 to 1.5 feet in 
thickness, while the overall thickness is approximately 90-100 feet. 
The upper submember can be seen at the Highways 401-15 and Highway 401 - 
Montreal Road roadcuts (photos 6a. and b. ) . 

Member B. . . . Member B has been sub-divided into two submembers, the 
lower submember having a thickness of approximately 20 feet and the 
upper having a thickness of approximately 40 feet. The contacts between 
members A and B and between the lower and upper submembers of Member B 
are exposed in a roadcut at Highway 401 and Montreal Road intersection 
east of Kingston, (photos 7a. and b. ) . 
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Photo 6a, 




Photo 6b. 



Photo 6a. and b. 



The upper submember of Member A of the Gull River 
Formation, taken at Highway 401 roadcuts, near 
Kingston. 
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Lower B 



Upper A 



- Note calcite vugs 



Phota 7a. Contact between Members A and B of the Gull River Formation 
at the Highway 401 and Montreal Road intersection. 




Photo 7b. Upper submember of Member B of the Gull River Formation west 
of Highway 401 and Montreal Road intersection. 




Photo 8. Member C at the highways 401 and 38 intersection, note the 
alternating beds of shale with lithographic limestone. 




Photo 9a. Stromatolite zones at the highways 401 and 38 intersection. 



13 










Photo 9b. Stromatolite zones at the highways 401 and 38 intersection. 




Photo 10. Member D exposed south of the highways 401 and 38 intersection, 
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The lower submember consists of thin to medium beds of gray litho- 
graphic limestone. This member is used as a building stone in Southern 
Ontario because of its uniform texture and hardness. In places it is 
almost pelletoidal, while at other localities it is thinly laminated. 

The upper submember ranges from a fine, semi-crystalline, litho- 
graphic limestone to a clastic, fossiliferous, argillaceous limestone. 
Its appearance varies locally; however, it can be distinguished from the 
underlying beds by its fossil content and from the overlying Member C by 
its lack of stromatolite zones. 

Member C. . . . This member, described in detail by Greggs and Sargent 
(1971) , is approximately 20 feet thick. It is essentially an alternation 
of thin to thick beds of dark gray, lithographic limestone and shale 
with massive beds of grey, semi-crystalline and lithographic limestone, 
which are locally fossiliferous (Photo 8) . This member is charcterized 
by four clastic carbonate zones containing stromatolites, the uppermost 
zone of which is continuous throughout the Kingston area (Peterson, 
1969) . The abundant fossil species were found by Greggs and Sargent to 
be time transgressive , thus difficult to use for correlation purposes. 
Structural features include mudcracks, cross-beds and ripple marks. 
Stromatolites beds and ripple marks are shown in photos 9a and 9b. 

Member D. . . . This uppermost member has been termed a "mixed" unit 
by Liberty (pers. comm.) on the basis of the wide range of grain sizes 
present in any one hand specimen (from sublithographic to finely granular 
limestone) . This member is composed essentially of massive beds (2-5 
feet), although thinner shaly beds are common. It is highly fossiliferous, 
locally containing yellow-orange flakes of oxidized material and weathering 
to a gray-buff colour, with a gritty texture. Member D, as well as 
Member C, is exposed in a roadcut at the highways 401-38 intersection 
(photo 10) . Although only about 15 feet of Member D is present in the 
basin, Liberty (1971) has recorded thicknesses of up to 55 feet in a 
roadcut east of Camden. The uppermost member has only been found in the 
Kingston and northern New York State areas and has not been traceable 
west of Napanee. 

Bobcaygeon Formation 

The Bobcaygeon Formation, as originally named and defined by Liberty 
(1963), is a lithogenetic unit lying between the lithographic limestone 
of the Gull River Formation and the interbedded limestones and shales of 
the Verulam Formation. The unit has been sub-divided into two members, 
15-feet and 5-10 feet thick, respectively. The lower member, exposed 
near the Camden Road - Highway 401 intersection, consists of a fine 
calcarenitic limestone with occasional beds of sublithographic limestone 
(Photo 11). The member is highly fossiliferous, the dominant fauna 
include Receptaculites occidentalis , Tetradium cellulosum , Dalmanella rogata , 
and Maclurites logani . 

The upper member, although equally fossiliferous, is less resistant 
to weathering and is often thick-bedded. It is essentially a gray, 
sublithographic limestone and claystone with a blue-gray, irregularly 
pitted, weathering surface. 

The Bobcaygeon Formation is biostratigraphically correlated with 
the Chaumont-Rockland-Leray-Hull strata, while lithologically the 
Formation can be traced into the uppermost Bony Falls Formation in 
Michigan and the Leray-Rockland beds of the Ottawa Formation in the 
Ottawa Valley (Liberty, 1969) . 
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Photo 11. Contact between the Gull River and Bobcaygeon formations 

as represented by beds of fossiliferous limstone overlying 
beds of sublithographic limestone. Photo taken near inter- 
section of highways 401 and 133. 




Photo 12. Verulam Formation exposed in an abandoned quarry in the 
area south of Morven. 



Verulam Formation 

The Verulam Formation, after Liberty (1955), is the youngest 
bedrock formation found in the Wilton Creek area. It is characterized 
by alternating bands of fine to medium-grained crystalline limestone and 
weathered, calcareous claystone, with shale partings (Photo 12) . Only 
the lower-most unit is exposed in the area. This unit is highly fossili- 
ferous with an abundance of "Prasoporid" bryozoans. The abundant fauna 
include: "Dalmanella" , "Prasopora" and echinoderms. Many of the fossils 
have been silicified, as can be seen at three abandoned quarries south 
of Chambers. The lack of exposure limits the control on bedrock contacts; 
therefore, the bedrock map (Map 1) in the southern part of the basin is 
relatively unchanged from Liberty (1971). Liberty (1971) estimated the 
thickness of the Formation to be about 350 feet in the Kingston area. 
The Verulam Formation is correlated with the Sherman Fall and lower-most 
Cobourg formations of New York State and the Ottawa Valley (Liberty, 
1971) . 

BEDROCK STRUCTURE 

The linear nature of the drainage basins and various structural 
features in southeastern Ontario appear to parallel the general trend of 
the Grenville Structural Province. The Precambrian rocks of the Grenville 
Province, have a northeasterly strike. Wynn-Edwards (1967b) suggests 
that the metasedimentary rocks in the Grenville have undergone a complex 
stage of folding and refolding resulting in a system of northeast- 
trending folds. These folds were subsequently reduced by erosion to a 
peneplained condition, the anticlinal fold axes remnant as shallow 
linear ridges. These ridges can be traced under the Paleozoic cover as 
is evident from scattered inliers. Major drainage features in the 
Precambrian can be traced upward into the Paleozoic rocks. Liberty 
(1971) suggests that the Precambrian surface topography is of major 
importance in the relief and distribution of the younger bedrock formations 
in the area. 

The Precambrian surface dips to the south, while the overlying 
Paleozoic formations dip south to southwest. The measured regional dip 
for the Simcoe Group in the Wilton Creek area is between 20 and 50 feet 
per mile, becoming steeper to the north. 

Several northeast- trending faults occur in the Napanee-Kingston 
area, as the result of post-depositional (Paleozoic) basement movements. 
The largest displacement (100 feet) has been recorded in the Napanee 
area (Liberty, 1971). Figure 3, a cross-section of the basin near Violet, 
presents the gamma-ray logs obtained from six MOE observation wells. 
From the logs, it is evident that the units are easily traceable across 
the valley, with two, 20-foot vertical displacements occurring, one on 
the north side and one on the south side of the creek. The northeasterly 
trend of Wilton Creek, in addition to these displacements, may suggest 
fault control of the creek. 
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The gamma-ray logs measure the inherent natural radioactivities of 
the rock units, with shales and shaly limestones generally showing a 
higher radioactivity than clean limestones and sandstones; therefore, in 
a gamma-ray log, the peaks in the intensity will correspond to shales, 
while the lows will indicate limestones, sandstones, etc. Digesting 
limestone samples correlating with high and low gamma values, in both 
dilute HC1 acid and ethylenediaminetetraacetic acid (EDTA) , indicate 
that the high radioactive samples have twice the shale (undigested 
remains) as the low samples. No attempt has been made to relate the 
gamma-ray logs to Liberty's stratigraphic sequence. At the present no 
continuous gamma-ray logs exist from the Precambrian upward into the 
Ordovician units in this area. 

Feature and joint orientation measurements for the Wilton Creek 
area are presented on the bedrock map (Map 1) . These orientations which 
were obtained in quarries and other outcrop exposures, are grouped 
together by area. There appears to be some correlation between the 
major fracture trends and the trend of Wilton Creek. Both the fracture 
trend and the trend of Wilton Creek are roughly oriented between 30° and 
50°. This suggests that the faulting was along fracture planes. 
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OVERBURDEN GEOLOGY 

The Wilton Creek drainage basin is blanketed by deposits of glacial, 
glaciof luvial and glaciolacustrine origin, representing deposition 
during the Wisconsinan stage of the Pleistocene Epoch. The glacial 
features, both erosional and depositional , indicate that ice movement 
was south-southwest. 

During deglaciation, meltwaters created a series of proglacial 
lakes in the Lake Ontario basin. These lake stages resulted in the 
formation of a number of 'abandoned beaches', marking stands in water 
levels. The best known stage, that of Lake Iroquois, is not represented 
in the Kingston -Napanee area. It is assumed that Lake Iroquois was 
ponded against the ice margin in this area. The lack of any substantial 
thicknesses of deposits, and the repetitive nature of some of the sediments 
in the Wilton Creek drainage basin, suggests deposition along a fluctuating 
ice front followed by a rapid lowering of the proglacial lake stages to 
the Lake Ontario stage. Local deposits of alluvium and muck of Recent 
age can also be found within the basin. The areal distribution of the 
Pleistocene and Recent deposits is indicated on Map 2. 

The overburden thicknesses, as shown on the isopach thickness map 
(Map 3) , averages about 3-5 feet on the bedrock highs and about 20-25 
feet in the bedrock valley. Local accumulations of sands and gravels in 
the bedrock valley have a maximum recorded thickness on the order of 90 
feet. 

A generalized stratigraphic sequence compiled from various roadcuts 
and sand pits within the drainage basin is listed below (Table 4) . 



TABLE 4. GENERALIZED STRATIGRAPHIC SEQUENCE IN THE WILTON CREEK AREA 



Stage 
Recent: 



Event or Formation 



Alluvium, marsh deposits 



Lithology 

gravel, sand, silt, 
clay, muck, peat 



Late 
Wisconsinan: 



Lake Belleville beach deposits 



gravel and sand 



Undifferentiated glacial, 
glaciof luvial and glacio- 
lacustrine deposits 



gravel, sand, silt, 
clay, till 



Upper 
Ordovician: 



Glaciof luvial deposits 
Simcoe Group 



gravel and sand 



limestone and shale 
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DIRECTION OF ICE MOVEMENT 

A number of glacial features can be used to determine the approximate 
direction of ice movement. In the Wilton Creek drainage basin, these 
features include striations, grooves and drumlins (Photo 13). The 
observed orientations are presented on the accompanying overburden 
geology map (Map 2) and are summarized below on Table 5. 



TABLE 5. GLACIAL FEATURE ORIENTATIONS 



Location 



Feature 



Trends 



Harrowsmith and north 
Harrowsmith - Wilton 



Striations 

Striations 
Drumlins 



234° - 242° 



226 

225° - 233° 



Wilton - Violet 



Striations 

Grooves 

Drumlins 



234° - 240° 

227° 

o o 
233 - 237 



Violet to lower 

basin boundary 



Striations 
Drumlins 



234° - 256° 
232° - 242° 



A visual comparison of the orientation of the creek valley resulting 
from faulting with the glacial orientations would appear to indicate 
that the preglacial topography added considerable control to glacial 
flow. Local striation directions strongly reflect the local valley 
trends (Map 2) . 




Photo 13. Glacial whaleback forms on Highway 401 near the Highway 133 
intersection. 






GLACIAL DEPOSITS 

The term "glacial deposit" is used here to refer to unconsolidated 
sediments deposited by direct glacial action and includes both till and 
stratified material incorporated within the till. The glacial deposits 
in the Wilton Creek drainage basin are thin and discontinuous, averaging 
about 3-5 feet thick in the north. 

The till, based on mechanical analyses presented later in the 
report, is typically sandy to silty with a higher clay content in the 
deposits in the southern reaches of the basin. It appears that the clay 
content in this area is from the incorporation of glaciolacustrine 
sediments into the till by the overriding of these sediments by the ice 
front. Three thin but stratigraphically distinctive till layers are 
present in the Wilton-Violet area. These layers are intcrbedded with 
glaciof luvial and glaciolacustrine sediments as is characteristic of 
deposition along an oscillating ice front. Similarly, clay balls can be 
found incorporated in the till near the contact between the till and 
glaciolacustrine clay. The interbedded Pleistocene deposits near Wilton 
are shown in Photo 14, while descriptive logs for roadcuts and pit 
exposures are presented in Appendix 1. 

Drumlins are present in the south and central portions of the 
basin. These drumlins, which belong to the Quinte drumlin field, are 
characterized by their low relief and elongated ridge-like nature (Mirynech, 
1967) . Photo 15 shows an elongated drumlin south of Morven. The average 
dimensions for length, width and height of the drumlins are 1000 x 300 x 
25 feet, respectively, while the average trend is 235°. 

Two basic types of drumlins are present; those composed entirely of 
till and those with a stratified sand and gravel core. The first type, 
being more abundant, can be found throughout the Napanee - Kingston 
area. The latter type has been described in detail by Deane (1950) and 
examples can be found locally within the drainage basin. A typical 
example of this type of drumlin is present southeast of Thorpe (Photo 
16). At this locality, a gravel pit is cut into the side of the drumlin 
exposing sand and gravel and about 3 feet of overlying till. Water well 
#1218, which has been drilled into this drumlin, encounters 10 feet of 
loam and 26 feet of sand and gravel. 

GLACIOFLUVIAL DEPOSITS 

The term "glaciof luvial" is applied to deposits of unconsolidated 
sediments (coarse silt to gravel) laid down by meltwater streams. These 
deposits are characterized by the well sorted and stratified nature of 
the materials. Both cross and graded bedding are common structural 
features related to these deposits. 

In the study area, extensive glaciof luvial deposits can be found on 
bedrock terraces in the creek valley. Photo 17 shows sand and gravel 
deposits on bedrock terraces near the community of Wilton. 

A subglacial origin is postulated for these glaciof luvial materials, 
the sand and gravel being deposited in ice caverns adjacent to and on 
bedrock terraces by subglacial meltwater streams. Thin layers of till 
and clay can occasionally be found interbedded with the sand and gravel. 
The largest sand and gravel deposit underlies Morven. Water-well logs 
indicate the presence of up to 90 feet of sand and gravel in a bedrock 
depression. Other large deposits of this type underlie the lacustrine 
clay near the mouth of Wilton Creek. 
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Sand 

Clay 

Till 
Sand 

Till 



Sand 



Photo 14. Interbedded till, clay and sand near Violet at the sanitary 
landfill site. 




Photo 15. Elongated drumlin south of Morven, 
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Photo 16. Drumlin with a sand and gravel core, in the Thorpe area. 
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Photo 17. Sand and gravel deposits on bedrock terraces near Wilton. 
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The glaciof luvial deposits near Harrowsmith and Wilton are presently 
being used as a source of sand and gravel for local road and building 
construction. Scott (1971) , in his examination of engineering properties 
of these deposits, obtained an approximate volume measurement of 100,000 
cubic yards for the one near Harrowsmith. He estimated the area of 
deposit near Wilton to be on the order of 1,000,000 square feet (Photo 
18) . Map 4 indicates the approximate and assumed extent of known sand 
and gravel deposits in the drainage basin. 

In the northern portion of the basin, two miles east of Hartington, 
a sand, gravel and till plug acts as the drainage divide between Wilton 
Creek and Knowlton Lake. This plug measuring about 600 feet by 200 feet 
completely fills the local bedrock valley (Figure 4) . A seismic refraction 
and electrical resistivity survey was carried out by the Ministry of 
Environment in order to determine the thickness of the overburden. The 
following stratigraphic sequence was interpreted from the results of 
this survey: 



Unit 



Depth Range (maximum) 
Resistivity Seismic 



sand or sandy till (dry) 
sand and gravel (wet) 
till or fractured bedrock 



0-5 feet 0-43 feet 

5-25 feet 

25 feet 43 feet 



The approximate thicknesses for the overburden at various stations 
on the plug are indicated in Figure 4. The maximum assumed thickness is 
43 feet. This plug was probably deposited during the final recession of 
the ice sheet. 

GLACIOLACUSTRINE DEPOSITS 



Glaciolacustrine deposits are defined as those deposits which were 
laid down in lakes adjacent to glacial ice masses, the materials being 
delivered to the lakes by meltwater streams. Particle sizes range from 
clay (deep water) to sand and gravel (nearshore and beach environments) . 

The glaciolacustrine clay and silt found interbedded with the 
glacial and glaciof luvial sediments in the Wilton Creek valley and 
adjacent low-lying areas are attributed to deposition in Lake Iroquois 
along a floating ice front. The clay in general become increasingly 
siltier towards the north as would be expected in a sedimentary progression 
(particle size increasing with decreasing distance from the glacial 
source) . The glaciolacustrine deposits vary in thickness from 2 to 5 
feet, the thicker beds occurring in topographic depressions. 

The clay and silt in the south-central portion of the valley can 
be subdivided into an upper and lower unit; the lower unit being more 
compact and darker in colour then the upper unit. Compression tests 
conducted on two similar clays in in the Kingston area suggest that the 
lower clay was glacially overridden (Sangrey, 1970) . This would support 
the concept of a fluctuating ice sheet in this portion of southeastern 
Ontario. 

Varved clay is only exposed in one area, the sanitary landfill site 
south of Violet. At this locality, glacial till is overlain by about 3 
feet of varved clay. The varved clay is overlain by another till layer, 
suggesting that its occurrence is the result of local ponding. Clay 
balls can be found incorporated in the upper till, near the till-clay 
contact. This may indicate turbulent water conditions prior to the 
advance of the ice front and subsequent till deposition. 
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Figure 4. Sand, gravel and till plug, east of Hartington. 
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Photo 18. Large sand and gravel pit near Wilton. 




Figure 6. Triangular plot of till samples collected in the basin. 
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The presence of "beach remnants" on several drumlins in the south- 
central portion of the basin has been attributed to a brief water stand, 
the Lake Belleville Stage (Mirynech, 1962). A pit cut into a gravel bar 
running off one of these drumlins (1*5 miles southeast of Morven) , 
indicates the presence of very coarse sand and gravel, with poorly 
developed bedding structure. The pebbles at this locality are well 
rounded and are primarily limestone and Precambrian fragments. The 
approximate elevations for these beach remnants are 375-425 feet. 

Scott (1971) suggests that the sand and gravel deposits in the 
creek valley near Wilton (400-foot elevation) are also beach remnants. 
He bases his suggestion on correlation of these deposits with other 
deposits in the Sydenham and Odessa areas which occur at the same 
elevation. Scott also found a number of clay balls in the sand and silt 
lenses of the Wilton deposits and attributed these to turbulent conditions. 
The writer feels however that this is not the case and instead postulates 
a subglacial origin for these deposits as is discussed previously. 

RECENT DEPOSITS 



The Recent deposits consist entirely of alluvium and muck. The 
alluvium deposits composed of silt and fine sand are very limited in 
area, lying adjacent to the creek. This silt and fine sand represent 
sediments eroded from upstream areas which have been reworked and sorted 
before being deposited along the stream bank. 

The muck or swamp deposits are much more extensive, comprising 
about 6% of the total basin area (MOE, unpublished land-use survey) . 
These deposits are confined essentially to topographic depressions, 
their poor drainage relating to both the impermeable underlying materials 
and their position relative to the local water table. These swampy 
areas are normally overgrown by water- tolerant plants, which upon decay 
add to the volume of the muck. 

The occurrence of swampy areas in the valley north of Harrowsmith 
is in part attributed to the high water-table conditions in the glacio- 
fluvial deposits flanking the valley sides and in part to the low perm- 
eability of the underlying Precambrian bedrock. The sands of these 
deposits are normally saturated year round, contributing to the maint- 
enance of the swamps by seepage. Other swampy areas in the southern 
part of the basin near the mouth of the creek are saturated alluvium 
deposits which have been subsequently overgrown by vegetation. The 
areal extent of these Recent deposits is presented on Map 2. 

MECHANICAL AND CARBONATE ANALYSES OF OVERBURDEN MATERIALS 

Mechanical and carbonate analyses were carried out on 15 overburden 
samples collected within the Wilton Creek drainage basin (Map 2). In 
addition, mechanical analyses were carried out on six sand and gravel 
samples collected from observation wells at the Violet landfill site. 
The mechanical analyses were evaluated using the Wentworth classifi- 
cation. 

Table 6 lists the results of these analyses along with a brief 
field description and the Trask (1932) sorting coefficient for each 
sample. Trask, considered a well sorted sediment to have a value of 
under 2.5. The mechanical analyses are plotted on a triangular diagram 
(Figure 5) . This graphic plot illustrates the limited grain-size 
composition of the till. 

Most of the samples collected south of Wilton appear to have a 
higher total carbonate and clay content but as the sample size was 
rather small no conclusions can be drawn. 
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TABLE 6. MECHANICAL AND CARBONATE ANALYSES OF OVERBURDEN MATERIAL 



Sample 

Number Field Description Gravel 

Samples Collected North of Wilton 
(sample #8 included for comparison) 



Grain Size % Sorting 
Sand Silt Clay (Trask) 



Calc. 



Dol. 



Total 
Carbonate % 



Ratio 
Calc. /Pol. 



(3) Till - sandy, reddish-brown 

(8) Till - silty, few pebbles, brown 

(9) Glaciof luvial sand - well sorted 

(10) sand and gravel 

(11) " sand and gravel 

(12) Till - sandy, silty, few pebbles 
(15) Till - sandy, brown 

Samples Collected South of Wilton 

(1) Till - sandy, many pebbles, compact, grey 

(2) Till - sandy, silty, brown 

(4) Till - many pebbles, compact, gray 
(upper till layer) 

(5) Till - many pebbles, compact, reddish- 
brown (middle till layer) 

(6) Till - sandy pebbles, compact, brown 
(lower till layer) 

(7) Till - sandy, silty, brown 

(13) Till - many pebbles, compact, brown 

(14) Till - sandy, brown 



30 
36 

10 



50 
34 
92 
70 
64 
50 
44 



34 
50 



40 



8 


2.8 


1.2 


2.8 


4.0 


0.4 


10 


3.2 


1.2 


7.0 


8.2 


0.2 


- 


1.7 


4.4 


45.8 


50.2 


0.1 


- 


2.1 


8.2 


39.8 


48.0 


0.2 


- 


2.1 


22.6 


18.6 


41.4 


1.2 


4 


2.0 


0.8 


1.6 


2.4 


0.5 


6 


2.8 


1.2 


2.2 


3.4 


0.5 



3 


25 


53 


19 


6 


42 


45 


7 


11 


29 


35 


25 


13 


33 


39 


15 


9 


39 


36 


16 


6 


37 


47 


10 


12 


48 


31 


9 


- 


58 


36 


6 



3.2 


10.8 


20.0 


30.0 


2.4 


21.0 


18.0 


39.3 


5.7 


22.0 


19.0 


41.0 


5.5 


15.2 


17.6 


32.8 


4.2 


18.2 


17.2 


35.4 


2.6 


15.6 


20.2 


35.8 


4.2 


1.0 


1.5 


2.5 


2.5 


21.0 


13.6 


34.6 



1.1 

1.2 



0.7 
1.5 



GLACIAL HISTORY 

No overburden deposits older than the late Wisconsinan are known in 
the Kingston-Napanee area. The glacial history as interpreted from the 
stratigraphic relationships observable in the study area is summarized 
below. 

The oldest deposits are of glaciof luvial origin, being deposited in 
the valley proper and along bedrock terraces flanking the valley walls. 
The presence of striae and grooves on the bedrock surface would suggest 
that earlier glacial sediments may have been present. The origin of the 
bedrock terraces is related to either or both glacial scour and stream 
processes. 

A subglacial origin is postulated for the glaciof luvial sediments, 
deposition occurring by meltwater streams flowing along the terraces and 
valley floor in ice caverns. Any earlier glacial sediments are assumed 
to have been eroded and incorporated into the stratified sand and gravel. 

With deglaciation, the waters of proglacial Lake Iroquois stood 
along the ice margin (Henderson, 1967) . Henderson, because of the lack 
of Lake Iroquois beach and nearshore deposits, suggests that the ice 
margin may have been floated clear of the land surface by the buoyant 
action of the Lake Iroquois waters. The type of the deposit which would 
be formed would therefore be dependent on the position of the ice front 
relative to where deposition was occurring. The presence of interbedded 
glacial, glaciof luvial and glaciolacustrine sediments in the south and 
central portion of the valley would suggest an oscillating ice front. 

The lower, glaciolacustrine clay layer in this area shows signs of 
being compressed, possibly by subsequent overriding by the ice front. 
In addition in the Violet and Wilton areas, clay balls have been found 
incorporated in the uppermost till near its contact with glaciolacustrine 
sediments. 

It is assumed that the ice sheet retreated rapidly from the area 
and was followed by a lowering of the lake level with the opening of the 
St. Lawrence outlet. The single lacustrine shore feature, a gravel bar 
south of Morven, has been related to a brief transitory post-Iroquois 
Lake stage, Lake Belleville (Mirynech, 1962) . 

The sand and gravel plug at the head of the creek valley can be 
attributed to the final stages of the glacial retreat from the basin. 



29 






GROUND WATER RESOURCES 

Ground water is the major source of domestic water supply in the 
Wilton Creek drainage basin. As such it is a valuable resource providing 
an inexpensive and near constant source of potable water. It is the 
intent of this report to evaluate this resource in terms of quantity and 
quality. 

INSTRUMENTATION 

A number of observation wells, streamflow gauging stations and 
meteorological stations have been established in the basin to provide 
geological, hydrological and climatological information. The instru- 
mentation present in the basin and its location are shown in Map 6. 
Hydrologic and climatologic data have been collected on a continuous 
basis since 1967. 

BEDROCK AND ITS HYDROGEOLOGIC CHARACTERISTICS 

Physical Characteristics 

The carbonate bedrock underlying most of the basin constitutes the 
major aquifer. Some basic assumptions inherent in the calculation of 
water-yielding properties of this aquifer, as based on its physical 
characteristics, are as follows: 

Ground-water availability and quantity . . . . Ground-water avail- 
ability and quantity is related to the number of interstitial openings 
in the bedrock and their degree of interconnection. The primary porosity 
(the natural volume of void space) of the carbonate rocks of the Simcoe 
Group is negligible as compared to their secondary porosity (fractures, 
solution channels and bedding planes). These openings, which act as 
interconnected conduits for ground-water flow, have been enlarged 
through the dissolution of CaC0 3 by undersaturated ground water (photos 

19 and 20) . Ground-water seepage from roadcuts and quarry faces in the 
area is confined to the above-mentioned openings and fractures. Photos 

20 and 21 show winter ice accumulations as the result of ground-water 
seepage from bedding planes. 

Well productivity The productivity of any specific well is in 

part dependent on the number of interstitial openings encountered, the 
amount of water available in the interstitial openings and the ability 
of these openings to transmit water to the well. Meinzer (1923), Lattman 
and Parizek (1964) and Siddiqui and Parizek (1969, 1971) concluded that 
wells situated on fracture traces in limestone bedrock, are considerably 
more productive than wells not situated in such positions. Similarly, 
wells in bedrock valleys are more productive then wells in bedrock 
upland areas. The observed variability of the productivity of wells 
drilled into the bedrock in different areas of the Wilton Creek basin is 
assumed to be related in part to variabilities in the density of fractures 
over the basin. 

Well productivity with respect to well depth Summers (1972) 

and Landers and Turk (1973) found that the number and size of secondary 
openings in crystalline bedrock and the associated productivity, decrease 
with increasing bedrock depth. Erosional unloading, in part responsible 
for surface fractures, and weathering which enlarges these fractures, 
has a decreasing influence on the bedrock with depth. A similar exami- 
nation of bedrock wells in the Wilton Creek basin supports the assumption 
of decreasing productivity with depth. 




Photo 19 




Photo 20 



Photos 19 and 20. Solution-widened fractures in the Wilton Creek 
drainage basin. 
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Photo 21 




Photo 22 



Photos 21 and 22. Ice accumulations as the result of seepage along 
bedding planes in the Wilton Creek area. 



Quantitative Characteristics 

The hydrogeologic properties of the bedrock can be expressed 
quantitatively by the coefficients of permeability, transmissibility and 
storage. These values were estimated for specific areas by pumping and 
recovery test data. 

Pumping and recovery tests . . . . These tests, in which the effects 
of pumping are measured in pumped wells and in nearby observation wells, 
were used here in arriving at reasonable estimates of the hydraulic 
properties of the bedrock. The data base consists of the results of 
three controlled tests conducted in the basin and the results of a 
series of pumping tests conducted at the Lennox Power Generating Station 
site southeast of the study area. This information has been supplemented 
by specific-capacity values obtained from short-term pumping tests which 
are discussed later. 

The Hantush and Jacob (1955), leaky artesian formula, as described 
by Walton (1960), was used in analyzing the pumping test data where 
observation-well measurements were available. The limitations imposed 
on this method by fracture porosity are discussed in Zeizel et al (1962). 
The Cooper and Jacob (1946), modified non-equilibrium formula, was used 
where observation-well measurements were not available. Table 7 presents 
the results of the pumping tests. The pumping-test data along with the 
method of computation are presented in Appendix II. The calculated 
coefficients of transmissibility range from 90-3900 Igpd/ft. A single 
permeability value is in the order of 100 Igpd/ft 2 and the coefficient 
of storage is estimated to be 4 x 10 - 4. 

In order to obtain a sufficient water yield to supply the needs of 
Service Centre E-5 located near the Highway 401 and Highway 144 intersection, 
it was necessary for the contractor to drill four separate wells. These 
wells were emplaced close together (radius of 3 feet) in an attempt to 
enhance fracturing of the limestone bedrock and thus effectively increasing 
the surrounding permeability. In addition, the last of these four 
wells, which is the producing well, has a larger diameter (8 inches) 
then the other three wells. Information on these wells is presented in 
Table 7. In the pumping test for the producing well (Figure 19) the 
fitted line was positioned in order to be compatible with the short-term 
character of the other tests and allow for comparison with these tests. 
It should be mentioned that the data indicate that a recharge boundary 
is reached after about 50 minutes and that a line fitted to the data 
points after this time period would result in a smaller drawdown and 
subsequently a higher transmissibility. 

The results of a series of pumping tests conducted at the Lennox 
Power Generating Station site southeast of the study area are included 
for comparison. The bedrock units are similar to those encountered in 
the study area. It was concluded in this study (Morrison, 1971) that 
the transmissibility varied from well to well, depending on the degree 
of interconnection of the rock fractures encountered. 

S pecific capacities .... Specific-capacity values calculated from 
the results of short-term production tests conducted on domestic wells 
by well drillers, were used to provide information on the hydraulic 
properties of the bedrock aquifer. Such data are not accurate or reliable 
and are used in this study only to obtain a rough approximation of 
aquifer characteristics. 
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TABLE 7. PUMPING-TEST RESULTS FOR WELLS FINISHED IN THE BEDROCK AQUIFER 



PUMPING RATE 
TEST WELLS (WELL I. D.) 



METHOD OF 
CALCULATION 



COEFFICIENTS OF 
PERMEABILITY TRANSMISSIBILITY STORAGE 



WC-7 
WC-4,WC-7 



3 Igpm (6") Hantush & Jacob (1955) -100 Igpd/ft 2 -1400 Igpd/ft 4 x 10 4 



3 Igpm (6") Cooper & Jacob (1946] 



Service Centre 

E5 100 Igpm (8"» Cooper & Jacob (1946,' 

Lennox Generating 

Plant 8 Igpm (6") NA 



»9S Igpd/ft NA 
-3900 Igpd/ft NA 
700-1400 Igpd/ft NA 



NA - Not Available 



Specific-capacity values are expressed as the yield of a well in 
gallons per minute per foot of drawdown (Igpm/f t) . To obtain a uniform 
data base these values were adjusted to represent a common well radius 
of six inches and a pumping interval of one hour. In order to reduce 
the error involved due to losses from storage, the data were selected to 
ensure that the amount of water withdrawn from the well greatly exceeded 
that in storage in the well. This required that values for wells which 
went dry during pumping be discarded. No adjustments were made to 
compensate for well loss due to turbulent flow in the vicinity of the 
well. 

The specific-capacity values can be used to obtain rough estimates 
of the transmissibility of the aquifer. Figure 6 illustrates the graphical 
relation between specific capacity and the coefficient of transmissibility 
as derived from Theis et al (1954) . Generally a high specific-capacity 
value is an indication of a high coefficient of transmissibility. The 
limitations of this method are documented in Theis' s paper. The estimated 
coefficients of transmissibility derived by this method for specified 
well depths are presented on Table 8. 

The adjusted specific- capacity values were divided by the total 
uncased portion of the well in saturated bedrock, to determine the 
specific capacity per foot of penetration or producitivity. Productivity- 
frequency graphs were prepared from these data by the method attributed 
to Kimball by Zeizel et al (1962) . This method assumes that the sample 
is log normally distributed. These graphs were used to examine the 
relationship between well yields and the following controls: 

1. wells finished within specific depth intervals in the bedrock, 

2. wells finished in bedrock having different lithologies, 

3. wells finished in either bedrock upland areas or bedrock valleys. 
The probable range of productivity for the different areas and 

depth intervals within the basin can be estimated from these graphs. 
The slope of the lines fitted to the plotted observation points is a 
function of the range in productivity. The steeper the line the wider 
is the range. Similarily the higher the position of the line on the 
graph the greater the productivity. 

The productivity-frequency graphs for wells finished within the 0- 
50- foot, 50.1-100-foot and the >100.1-foot depth intervals are illustrated 
in figures 7, 8 and 9. In these graphs, the productivity data have been 
arbitrarily divided into three groups from northwest to southwest- the 
Township of Portland, the Township of Ernestown and the townships of 
North and South Fredricksburgh. These groupings roughly coincide with 
changes in the bedrock lithology at the surface, older units outcropping 
in the northeast and younger units outcropping in the southwest (Map 1) . 
Because the bedrock units dip to the southwest, the unit encountered 
within a specific depth interval in one township will not necessarily be 
the same as that encountered within the same depth interval in the other 
townships. Very few wells in the townships of North and South Fredericks- 
burgh were completed to a depth of >100.1 feet. The limited size of the 
well population precluded the preparation of a graph for this depth 
interval. 

In the Township of Portland, wells completed to a depth of >100.1 
feet are finished in either Precambrian or Paleozoic bedrock depending 
on the well location. 

In order to compare the productivities of wells finished in different 
bedrock units these wells are grouped by the unit in which they are 
finished. The productivity- frequency graphs illustrated in Figure 9 
present the results of this grouping. The productivity-frequency graphs 
illustrated in Figure 10 compare the productivity of wells finished in 
the bedrock valley and upland areas. 
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Figure 6. Relationship between specific capacity and the coefficient of transmissibility. 
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Figure 7. Productivity-frequency graphs for bedrock wells completed in the 
<60 foot depth interval. 
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Figure 8. Productivity-frequency grephs for bedrock wells completed in 
the 50.1-100 foot depth interval. 
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Figure 9. Productivity-frequency graphs for bedrock wells completed in the > 100.1 foot 
depth interval. 



39 









0.01 



0.0001 







\ 


1 
















LEGEND 

Bedrock vallev 

— •— • Badcock upland 








\ 


1 


































' 


\ 

\ 




\ 
































"v| 


1 

V 

1 


\ 

k 

i 


































1 




\ 
\ 

>- 


> 

\ 


\ 


■*. 
































> 


i 









0.01 0.06 0.2 0.5 1 2 5 10 20 30 4050 60 70 80 90 95 98 99 99.8 99.99 

0.1 99.5 99.9 



Percentage of wells with productivity equal to or more than stated value 

Figure 10. Productivity-frequency graphs for wells finished in bedrock 
valleys and uplands. 



TABLE 8. TRANSMISSIBILITIES ESTIMATED FROM SPECIFIC-CAPACITY DATA FOR WATER WELLS 



Range of 
Penetration (ft) 



Calculated Mean 
Specific Capacities 



Carbonate Bedrock 

0-25 
25.1-50 

50.1-75 
75.1-100 

100.1-100 

125.1-150 
150.1+ 



Average for Bedrock Wells 0.912 
Average for Overburden Wells 1. 18 



Average for Overburden and 
Bedrock Wells Combined 



Total Number 
of wells 



% of Total 
wells 



Estimated Transmissibility 
from Graph (Figure 6) 



1.39 


70 


12.96 


3,600 


1.14 


195 


36.11 


3,000 


0.43 


144 


26.67 


1,075 


1.41* 


65* 


12.04 


3,650 


0.46 


63 




1,150 


0.58** 


28** 


5.18 


1,450 


0.07 


27 




160 


0.075 


24 


4.44 


170 


0.084 


14 


2.59 


190 


0.912 


(Total) 540 


(Total) 99. 99 


2,300 


1.18 


27 




1,400 


0.924 


(Total) 567 




2,400 



Igpd/ft 



* Two anomalous samples (45.6 and 53.17 Igpm/ft) 
** One anomalous sample (14.28 Igpm/ft) 



The number, minimum, maximum, geometric mean and the median values 
of productivity for wells situated in different hydrogeological settings 
are presented in Table 9. These data were obtained from the productivity- 
frequency graphs and the statistical treatment of the productivity 
values within the categories. The following trends are apparent from an 
examination of the data presented in the graphs and on Table 9: 

1. The Productivity of wells finished in the same townships (same 
area) decreases with increasing well depth. In other words, 
shallow wells appear to be on the whole more productive then 
deep wells. 

2. The productivity of wells finished in the bedrock valley are 
greater than those in the bedrock uplands. A similar study by 
Siddiqui and Parizek (1971) , found this relationship to hold 
true for carbonate rocks in central Pennsylvania. They attri- 
buted the higher productivity in wells finished in the bedrock 
valleys to a greater concentration of fractures in the valleys. 
It is suggested here that another factor for the higher produc- 
tivity found in the Wilton Creek basin may result from the 
fact that because the water level is much closer to the ground 
surface in the valley than in upland areas, the completed 
depth of the wells finished in the valley would on a whole be 
less than that in the uplands area. 

3. The productivity variations of wells finished at set depth 
intervals in areas having different lithologies are not readily 
apparent. 

The Kruskal-Wallis nonparametric statistical test as described by 
Siddiqui and Parizek (1972) was applied to the productivity data in 
order to determine whether the observed differences in well productivities 
for areas having different bedrock lithologies are significant or random. 
The results of this test are presented on Table 10, the statistic used 
(H) is defined by the following formula: 

H = 12 E £ " 3(N+1) 
N (N + 1) j = l nj 

where: k = number of samples, 

nj = number of observations in the j th sample, 

N = Sn j , the total number of observations, 

Rj = sum of the ranks in the jth sample (column) , 
k 

I = sum over the k samples (columns), 
j = l 

For nj>5, as is the case, H is distributed approximately as chi 
square (x2) with degrees of freedom (df ) = (k -1 ) . 

The mull hypothesis to be tested is that there is no significant 
difference in the productivities of wells finished within areas of 
different bedrock lithologies. The alternate hypothesis is that the 
productivities are significantly different. Setting a significance 
limit of a = 0.05 with the degrees of freedom (df) - 2, the significance 
level for chi square corresponding to H is 5.991. All calculated values 
of H which are larger than 5.991 indicate that the observed differences 
in productivities are significant. Those values smaller then 5.991 
indicate that the differences are random or not significant. Based on 
this statistic (Table 10), the productivities of wells in the three 



TABLE 9- PRODUCTIVITY VALUES (IGPM/FT/FT) OF WELLS COMPLETED IN DIFFERENT HYDROGEOLOGICAL SETTINGS 

n minimum maximum median geometric mean 



Bedrock wells finished within specified depth intervals 



<50.0 Feet 
50.1 - 100.0 Feet 
>100.1 Feet 



272 



2.0 x 10 
2.9 x 10 
2.7 x 10 



-4 
-5 
-5 



8.0 

1.3 

7.1 x 10 



-2 



NA 



1.7 x 10 
2.9 x 10 
1.7 x 10 



.-1 
-2 
-3 



Wells finished in specific rock types 

Verulam Fm. - Bobcaygeon Fm. (N.S. Fredricksburg, <50.0') 
Bobcaygeon Fm. - Gull River Fm. , Mbrs.D,C,B 

(N.S. Fredricksburgh, 50.1-100.0') 

Gull River Fm. , Mbrs. D,C,B upper (Ernestown, <50.0') 

Gull River Fm. , Mbr. B lower, Mbr. A upper & middle 

(Ernestown, 50.1-100.0') 
Gull River Fm., Mbr. B lower, Mbr. A upper & middle 

(Portland, <50.0') 

Gull River Fm. , Mbr. A middle & lower (Ernestown, >100.1') 

Gull River Fm. , Mbr. A middle & lower 

(Portland, 50.1-100.0') 

Gull River Fm., Mbr. A middle & lower (Portland, >100. I 1 ) 
Shadow Lake Fm.,- Potsdam Fm. - Precambrian 

(Portland, >100.1') 



31 

17 
110 

77 



3.0 x 10 



-4 



1.3 x 10 
3.0 x 10 



-4 
-2 



1.1 x 10 



-4 



131 


2.0 x 


10 


7 


9.0 x 


10 


118 


2.9 x 


10 


40 


2.7 x 


10 


19 


3.0 x 


10 



-4 
-5 



-5 

-5 



-5 



8.0 

1.9 x 10 
4.1 

9.6 x 10 

2.6 

2.7 x 10 

1.3 
4.1 x 10 

7.1 x 10 



-2 



-1 



-3 



-3 



-2 



2.7 x 10 



-2 



8.0 x 10 

1.1 x 10 



-4 
-2 



1.8 x 10 



-3 



1.9 x 10 
1.5 x 10 



-2 

-4 



2.1 x 10 
2.8 x 10 



-3 

-4 



3.1 x 10 



-4 



5.9 x 10 



-1 



2.3 x 10 
1.5 x 10 



-3 
-1 



4.1 x 10 



-2 



8.2 x 10 
5.6 x 10 



-2 
-4 



2.6 x 10 
5.3 x 10 



-2 
-3 



4.5 x 10 



-3 



Wells finished in bedrock uplands and bedrock valleys 

Bedrock Uplands 
*. Bedrock Valleys 



3.1 x 10 
1.4 x 10 



-5 
-4 



9.9 x 10 
2.6 



-1 



1.9 x 10 
5.5 x 10 



-3 
-3 



NA - not available 



4.3 x 10 
9.3 x 10 



-2 

-2 






' 



(r^ 






TABLE 10. KRUSKAL-WALLIS ONE-WAY ANALYSIS OF VARIANCE FOR THE SPECIFIED PARAMETERS 



Parameters Examined 



Degree of 
Freedom 


Calculated 
H 


2 


7.26 


2 


7.22 


2 


0.99 



Significance 
Level 2 x 0.05 



Variation 
Significance 



For wells grouped by Townships 

Finished in <50.0 foot range 

Finished in 50.1-100 foot range 

Finished in > 100.1 foot range 

Wells finished in the bedrock valley 
as compared to bedrock uplands 



8.55 



5.991 


Significant 


5.991 


Significant 


5.991 


Not Significar. 


3.841 


Significant 



townships are significantly different for the 0-50 foot and 50.1-100 
foot ranges of bedrock penetration. There is no significant difference 
in wells finished in the > 100.1 foot range even though the bedrock 
lithologies vary from Precambrian rocks, to the 'basal group' elastics 
to limestones. Therefore the null hypothesis is rejected for the wells 
finished in the 0-50 foot and 50.1-100 foot ranges and accepted for the 
> 100.1 foot range. 

A Kruskal-Wallis test, conducted on the productivities of the 
upland and valley wells, also indicates that there is significant 
difference in productivity* df=l, significance level 3.841 (Table 10). 

OVERBURDEN DEPOSITS AND THEIR HYDROGEOLOGIC CHARACTERISTICS 



The sparseness of overburden in the drainage basin is illustrated 
by the overburden (isopach) thickness map (Map 3) . The average thickness 
of the overburden, with the exception of local sand and gravel deposits 
in the Wilton Creek valley proper, is in the order of 3-5 feet (Photo 24) . 

The till and lacustrine deposits act as a low-permeable mantle, 
creating semi-confined aquifer conditions where they overlie glaciof luvial 
sediments and bedrock. The sands and gravels of the glaciof luvial 
deposits locally constitute favourable aquifers; however, their thickness 
and areal extent are too small to yield anything but small domestic 
water supplies. The areal distribution of these glaciof luvial deposits 
is presented in Map 4. From the water-well records, the maximum thickness 
of these deposits is about 90 feet, while the average thickness is about 
20 feet. Photo 25 shows de-watering operations in a gravel pit near 
Wilton. 

The mean specific-capacity value (1.18 Igpm/ft) for 27 overburden 
wells was used to estimate a mean coefficient of transmissibility of 
1,400 Igpd/ft, with a coefficient of storage in the order of 0.1 (Figure 11). 

Estimates of the permeabilities of various glaciof luvial sands 
exposed in sand pits in the Wilton bedrock valley were obtained using 
Hazen's empirical relationship between grain-size distribution and 
permeability (Hough, 1957) . The results are listed below: 






Sample 



11 

Well #5(10-15') 
Well #6(12-13') 
Well #7(6-8' ) 
Well (12*) 
Well (22-23') 
Well (28-30') 



Location 

East of Hartington 
East of Hartington 
East of Hartington 
Violet Landfill 
Violet Landfill 
Violet Landfill 
Violet Landfill 
Violet Landfill 
Violet Landfill 



10% Finer (mm) Permeability (Igpd/ft ) 



2.3 x 10~ 

3.5 x 10 

4.6 x 10~ 
1.3 x 10~j" 

5.5 x 10~ 
1.65 x 10~ 

2.6 x 10~ 

2.1 x 10 

1.2 x 10 



550 

850 

1000 

400 

3000 

450 

600 

510 

375 



Mean = 860 Igpd/ft' 




Photo 23. Sparse overburden in the Violet area. 




Photo 24. Dewatering of saturated sand and gravel in a pit near Wilton. 



GROUND WATER MOVEMENT 

Ground water moves continuously through the hydrogeological environ- 
ment. The rate of movement is dependent on the local hydraulic gradient 
and the hydraulic properties of the medium. The theory of ground-water 
movement within the three-dimensional domain has been described in 
detail by Hubbert (1940), Toth (1962, 1963) and Freeze and Witherspoon 
(1966, 1967). 

In theory, ground-water movement is confined on all sides by 
imaginary, impermeable boundaries. In the Wilton Creek basin, the 
lateral boundaries as established from water-well records appear to 
coincide with the drainage divide. The lower boundary underlying the 
domain is set at a horizon where flow approaches static conditions. 
This horizon is arbitrarily set at the elevation of the streambed near 
the mouth of Wilton Creek (=200 foot m.s.l.) and is extrapolated under 
the rest of the basin. This presupposes no underflow into nor out of 
the domain and that all ground-water discharge is into the stream. The 
upper boundary is represented by the water table, the upper surface of 
the zone of saturation. 

Hydraulic Potential 

The hydraulic potential of the ground-water at any point in the 
domain is defined by Hubbert (1940) as follows: 

$ = g z + (p - P Q )/P 

where: 

$ = hydraulic potential at any point in the domain, 

g = acceleration due to gravity, 

z = elevation at the given point above a standard horizontal 

datum, 
p = atmospheric pressure, p = pressure at the given 

point, 
p = density of water. 

The hydraulic head ($) or level to which the ground-water will rise 
in a well or piezometer (measured in feet of water above a standard 
datum) , equals the hydraulic potential divided by the constant (g) . The 
water-level map (Map 5) , which illustrates the relative values of hydraulic 
head as measured in various water wells in the basin, can be used to 
approximate the horizontal direction of ground-water movement (horizontal 
flow vector) . Ground-water movement is from areas of high hydraulic head 
on this map to areas of low head, in a direction perpendicular to the 
water-level contours. From the configuration of the water-level contours, 
it is evident that the local topographic relief has a strong influence 
on the hydraulic head, the water-level map being a subdued replica of 
the topographic base map. Ground-water movement is directed down the 
local topographic slope and in a more general pattern from the top end 
and divides of the basin to the mouth. 



The vertical direction of ground-water movement and the magnitude 
of the hydraulic gradient for a specific location can be established by 
measuring the hydraulic head in two or more piezometers which are placed 
side by side and set at different depths. In recharge areas the measured 
hydraulic head decreases as the piezometer depth increases. In discharges 
areas, the reverse holds true. Hydrographs for two piezometer nests, 
one (WC-1) located in the northern and the other (WC-8) in the southern 
reaches of the basin are presented in figures 11 and 12. 

Piezometer nest WC-1 has two piezometers; one set at 25 feet and 
the other set at 63 feet. Based on the measurements for 1968 to 1973, 
ground-water movement is vertically downward. This suggests that this 
portion of the basin acts as a recharge area. Piezometer nest WC-8 also 
has two piezometers; one set at 29 feet and the other at 176 feet. 
Based on the measurements, ground-water movement is vertically upward as 
is characteristic of a discharge area. 

Hydraulic Properties of the Medium 

The hydraulic properties of the overburden and bedrock aquifers 
have been discussed previously. A crude estimate for an average 
coefficient of transmissibility for the basin is obtainable by applying 
Darcy's Law as follows: 

Q = 2 TIL 
where: 

Q = mean baseflow discharge for the period under investigation (1967- 
74) which is approximated as 12.7 cfs, 

T = average transmissibility, 

I = average water - table slope perpendicular to the stream which is 
estimated as 0.02 feet/foot, 

L = creek length north of the recording station which is approximately 
21.5 miles. (Length is multiplied by 2 to simulate flow into the 
stream from both sides.) 

The coefficient of transmissibility thus calculated is 1360 Igpd/ft . 
This value compares favourably with values estimated from pumping tests 
(Table 7) and specific-capacity data (Table 8) . 

WATER LEVEL FLUCTUATIONS 

Water-level measurements have been collected at regular intervals 
from the eight MOE - IHD observation wells in the basin. The water 
levels, which fluctuate continually, reflect the dynamic nature of the 
ground-water storage in the basin. This fluctuation is related to 
changes in the ratio of ground-water recharge to discharge. A rise in 
the ground-water level implies that recharge exceeds discharge and water 
is taken into storage. A decline in the ground-water level implies that 
the reverse holds true. 
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Figure 11. Hydrograph for MOE observation well WC-1 and the precipitation record for 
Hartington for 1968 to 1973. 
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Figure 12. Hydrograph for MOE observation well WC-8 and the precipitation record for 
Morven for 1968 to 1973. 
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TABLE 11. SUMMARY OF WATER-LEVEL FLUCTUATIONS IN THE BASIN 







1967- 


68 


1968- 


-69 


1969 


-70 




1970- 


-71 


1971 


-72 




1972 


-73 




Piezometer Max. 




Max. 




Max. 






Max. 




Max. 






Max. 




Well 


or well 


and 




and 




and 






and 




and 






and 




Number 


Depth 


Min. 


Var. 


Min. 


Var. 


Min. 


Var. 




Min. 


Var. 


Min. 


Var. 




Min. 


Var. 


WC-1 A 


25 


5.5-15.7 


10.2 


4.9-15.8 


10.9 


4.5-16.0 


11.5 


3 


1-16.7 


13.6 


5.1-15.8 


10.7 


4 


7-15.0 


10.3 


WC-1 B 


63 


18.0-27.5 


9.5 


11.9-26.9 


14.9 


12.2-27.7 


15.5 


8 


6-28.2 


19.6 


16.7-28.4 


11.7 


13 


7-26.3 


12.6 


WC-2 


78 


28.1-55.6 


27.5 


26.1-51.8 


25.7 


28.6-57.5 


28.9 


23 


8-57.8 


34.0 


25.3-55.5 


30.1 


27 


6-57.5 


29.9 


WC-3 A 


178 


21.1-22.9 


1.8 


19.3-24.0 


4.7 


24.3-31.7 


7.4 


23 


4-29.7 


6.3 


25.5-29.1* 




23 


0-28.8* 




B 


118 


26.7-33.6 


6.9 


20.4-36.9 


16.5 


23.8-32.1 


8.3 


24 


1-32.7 


8.6 


28.9-32.7 


3.8 


25 


0-31.1 


6.1 


C 


46 


18.9-35.5 


6.6 


7.7-36.1 


28.4 


11.3-36.2 


25.9 


6 


1-37.1 


31.0 


16.8-35.5 


18.7 


7 


4-34.3 


26.9 


D 


29 


11.2-26.7 


15.5 


5.8-28.5 


22.7 


8.7-29.1 


20.4 


5 


7-33.0 


27.3 


12.9-33.9 


22.0 


6 


4-22.6 


16.2 


WC-4 


44 










2.3-9.5 


7.3 


4 


6-8.1* 




3.6-8.5* 




3 


3- 8.7 


5.4 


WC-5 


103 


4.1-20.4 


16.3 


3.4-20.6 


17.2 


2.9-20.9 


18.0 


3 


4-20.8 


17.4 


3.7-17.2 


13.5 


3 


5-20.5 


17.0 


WC-6 (1)C 


85 


33.6-37.0 


3.4 


29.6-37.7 


8.1 


29.0-36.7 


7.7 


29 


3-36.9 


7.6 


33.1-36.9 


3.8 


28 


9-35.6 


6.7 


(2)B 


45 


19.6-30.9 


11.3 


18.0-23.4 


5.4 


16.9-22.4 


5.5 


18 


0-23.6 


5.6 


20.8-24.3 


3.5 


18 


0-22.1 


4.1 


(3)A 


13 


10.1-12.0 


1.9 


9.9-12.9 


3.0 


10. 1-11.9 


1.8 


10 


3-11.8 


1.5 


10.4-12.0 


1.6 


9 


3-11.2 


1.9 


WC-7 


23 


13.0-19.0 


6.0 


12.3-18.8 


6.5 


12.1-18.5 


6.4 


8 


3-16.9 


8.6 


8.1-16.1 


8.0 


8 


0-15.1 


7.1 


WC-8 FCS 


176 


1.0- 2.6 


1.6 


0.1- 2.5 


2.4 


0.1- 3.1 


3.0 




F- 2.6 


2.6+ 


1.6- 2.3 


0.7 





3- 5.7 


5.4 


RDS 


29 


5.6- 7.1 


1.5 


2.4- 7.6 


5.2 


2.1- 7.9 


5.8 


1 


7- 8.3 


6.4 


2.95- 7.7 


4.75 





3-6.1 


5.8 


* Partial 


Record 






























F = Flowing 








































The fluctuations of the water level at two locations in the basin 
are represented graphically by hydrographs for the observation wells in 
figure 11 and 12. These hydrographs indicate that the water levels 
fluctuate seasonally, with lows occurring in the late summer and highes 
during the spring melt. Similarly both the magnitude of the fluctuations 
and the direction of flow, as represented by the hydraulic head difference 
in closely spaced piezometers, varies with location. A summary of the 
seasonal water-level fluctuations is presented in Table 11 for different 
observation wells in the basin. 

HYDROLOGIC BUDGET 

The hydrologic budget, a quantitative statement of the balance 
between recharge and discharge for a set interval of time, can be 
expressed mathematically by the following equation: 



P = R +R +ET + U + S 
g s - - 



where : 



P = total precipitation, 

Ra = ground-water runoff or baseflow 

Rs = surface runoff, h + R s = R t = total runoff 

ET = evapo transpiration, 

U = subsurface underflow, 

S = change in storage (surface, soil moisture, ground water) 

S = ( S. + S + S ) 

s sm g 

This equation is calculated by estimating and measuring the various 
parameters which affect its solution. Precipitation and total runoff 
are measured directly, the surface and ground-water runoff components 
and the evapotranspiration are estimated. It is assumed from the water- 
level map (Map 5) that the subsurface underflow into and out of the 
basin is negligable. Similarly for the period under investigation 
(1967-68 to 1973-74 water years) it is assumed that there is no appreci- 
able change in storage. The equation can therefore be simplified to: 



P = R + R + ET 
g s 



Precipitation 



Precipitation which includes rain, snow and dew is measured directly 
at the meteorological stations in and around the basin. Precipitation 
and temperature records for these stations were generally found to be 
erratic; however, two stations (Hartington IHD and Morven IHD) had 
sufficient data to cover the seven-year period. Mean monthly precipitation 
and temperature data tabulated for these two stations are presented in 
Table 12. 

Runoff 

Runoff is that portion of precipitation which appears in surface 
streams and is measurable by water-level gauges installed on the stream. 
Data from the hydrometric station 2HM-004, near Morven, were used to 
approximate the total stream runoff from the basin. 



The ground-water runoff component or baseflow, was estimated from 
the hydrometric records by the method attributed to Wundt by Meyboom 
(1967). Wundt' s method of analysis assumes that the monthly minimum 
streamflow is taken as a measure of baseflow. The annual baseflow for 
each water year (Table 12) is the mean monthly minimum flow for that 
year. This simplistic method provides a means of applying a statistical 
rather than an arbitrary approach towards baseflow analysis. 

Evapo transpiration 

It is necessary to measure cvapotranspiration by indirect methods. 
Two different methods were incorporated into this report, the Munson- 
Index Method (ASCE, 1966) and the Water Balance Equation or Indirect 
Method. 

The Munson-Index Method is an empirical method in which the potential 
evapotranspiration is taken directly from a table on which monthly mean 
temperatures are compared to potential evapotranspiration values. This 
method assumes that there is no evapotranspiration below the monthly 
mean temperature of 28.4°F. The empirical results are presented on 
Table 12. 

The Indirect Method consists of solving the simplified water 
balance equation: 

ET = P - P^. 

The results are presented on Table 12. 

In comparing these two methods, the results obtained by the Munson- 
Index Method are approximately 10% higher than those obtained by the 
Indirect Method. The difference in these two methods can readily be 
explained by the cumulative errors inherent in the estimation of the 
various components of the budget, and also the net changes in basin 
storage and underflow which were assumed to be negligible. The seven- 
year average values for both Munson-Index and the total calculated lake 
evaporation (Hartington-IHD station pan evaporation measurements) are 
also in close agreement (Table 12) . 

These methods do not take into consideration all of the factors 
which contribute to the actual evapotranspiration such as soil types, 
vegetation types, air humidity, wind velocity; however, they probably 
give a valid approximation. 
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TABLE 12. ESTIMATED ANNUAL HYDROLOGIC BUDGET OF THE WILTON CREEK DRAINAGE BASIN FOR THE INTERVAL 1967-68 AND 1973-74 



Water Year 



Precipitation* 
Inches 



Temperature F* 



1967-68 
1968-69 
1969-70 
1970-71 
1971-72 
1972-73 
1973-74 
Average (7 Years) 



37.34 E 
32.78 E 
33.96 
36.45 E 
42.86 
39.40 
36.45 
37.04 



43.67 
43.80 
42.94 
42.97 
43.67 
44.89 
43.57 



* Average for Hartington and Morven 

E Value based on some estimated daily values 

Potential Evapotranspiration 

Baseflow as % Munson P.E. Indirect 

of Precipitation Index Method Method 



10.1 
13.6 

9.1 
12.8 

8.5 
11.9 

9.9 
10.7 



23.24 
22.26 
22.81 
22.36 
22.13 
23.62 
22.33 
22.68 



20.82 
13.01 
20.21 
23.13 
28.43 
17.57 
21.25 
20.64 



Total Runoff 


Estimated 


Baseflow as % 


Inches 


Baseflow 


of 


Total Runoff 


16.52 


3.78 




22.9 


19.76 


4.45 




22.5 


13.76 


3.10 




22.5 


13.32 


4.67 




35.1 


14.44 


3.64 




25.2 


21.84 


4.68 




21.0 


15.20 


3.60 




23.7 


16.40 


3.97 




24.2 



Pan . Evap . 
(Hartington) 



18.04 E 
22.49 
19.77 E 
23.72 E 
23.27 
26.83 
23.78 
22.56 E 






HYDROCHEMISTRY 

The hydrochemistry is an important consideration in the evaluation 
of the water resources of an area. The water quality characteristics 
not only dictate the suitability of the waters for use but also can be 
useful in estimating the general ground-water flow patterns. 

As the waters move through the drainage basin they are subject to a 
series of continuous reactive processes which tend to change their 
chemical composition. These changes either contribute additional ions 
or remove them depending on the chemical reactions which are involved. 
Specific trends emerge in the hydrochemistry when viewed on either a 
regional or a basin-wide scale. For instance ground water in a recharge 
area generally has a bicarbonate character while that in a regional 
discharge area has a chloride character (Chebotarev, 1955) . In the 
following sections the hydrochemistry of the basin's waters will be 
presented and discussed with respect to their suitability for domestic 
use. In addition the hydrochemical data is used to estimate general 
ground-water flow patterns. 

Sampling Procedures and Analyses 

A total of 38 water samples were collected from wells in the Wilton 
Creek drainage basin, from October 16-19, 1973 to establish the chemical 
characteristics of the ground water. In addition, one sample collected 
on November 16, 1973 and 14 samples collected during a water survey 
conducted in 1969 are included in the presentation. Both MOE observation 
wells and private domestic wells were sampled. 

Sampling procedures varied with the type of well. The domestic 
wells (normally equipped with electric or hand pumps) were pumped for a 
five minute period of time in order to ensure fresh samples. No water 
samples were taken from wells where a chlorinator or water softener 
could not be by-passed. The observation wells and piezometers were 
sampled by bailing these dry where possible and then allowing them to 
recover prior to sampling. In a few cases, the well depth was excessive 
and the well could not be bailed dry. In such cases the well was bailed 
several times allowing the water level to recover between bailings 
before samples were taken. 

All water samples were analysed for bicarbonate alkalinity, pH, 
temperature and conductivity at the time of sampling. In addition, 
static water levels were recorded for the MOE wells. Laboratory analyses 
included: calcium, magnesium, sodium, potassium, chloride, iron, bicarbonate 
alkalinity, sulphate, nitrate, total phosphorus, total dissolved solids, 
hardness and alkalinity as CaC03, (all in milligrams/litre) and pH and 
conductivity (micromhos per cm) . The results of the water analyses are 
presented in Table 13 (in Appendix C) . 

Water samples had been taken periodically from the creek during the 
early stages of the IHD representative basin study, to establish background 
information on stream water chemistry. Table 13 present the results of 
these analyses for November 1968 to April 1969. Only minor variations 
in water quality occurred over that period. 
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A computer program written by Wigley (1971) and modified by J. 
Coward (pers. comm. ) was employed to obtain mean maximum and minimum 
relative (ion ratios) and absolute values for different chemical con- 
stituents in the water. These results, in addition to saturation values 
of calcite, dolomite, gypsum and magnesite, are included in Table 15 (in 
Appendix C) . 

The saturation indices for calcite (SATCAL) , dolomite (SATDOL) , and 
gypsum (SATGYP) , as defined by Wigley (1971), were determined as follows: 

SATCAL = log(aCa ++ aHCO~ K HCQ )-log(aH + K ), 

SATDOL = log(aCa ++ aM ++ (aHCC~) 2 (K HC0 ) 2 -log ( (aH + ) 2 K ), 

SATGYP = log(aCa ++ aSO~~)-log(K 



4 »~y 



where: 



++ 
aCa is the activity of the calcium in solution, 

++ . 
aMg is the activity of the magnesium in solution, 

+ 
aH is the activity of the hydrogen in solution, 

aHCO is the activity of the bicarbonate in solution, 

aSO~" is the activity of the sulfate in solution, 

K HCO-j' K c ' K ^' K „ anc * K „ are temperature dependent equilibrium constants. 

j ^ Q Cj III 

In addition Coward (pers. comm.) included the saturation index for 
magnesite (SATMAG) as follows: 

SATMAG = log(aM ++ aHCO" K.,__ )-log(aH + K ) 

The saturation indices were calculated using temperature dependent 
equilibrium constants given by Wigley and the lab values of temperature 
and pH. The calculated indices of saturation are considered to be 
significant only if the values deviate from zero (the saturation point) 
by 0.4 or more (the analytical error through using lab values, J. Coward, 
pers. comm.). Positive values (greater then 0.4) are considered to be 
supersaturated while negative values (less then -0.4) are considered to 
be undersaturated. 

The saturation indices are indicative of the ability of the water 
to dissolve the different mineral species such as calcite or dolomite 
from the aquifer materials. For example, where the ground water is 
undersaturated with respect to calcite the water has the potential for 
dissolving calcite. Conversely, where the water is saturated with 
calcite, calcite may be precipitated out if C0 2 is lost. 

The pumpage of a well will result in a drop in aquifer pressure in 
the vicinity of the well. If the ground water in the vicinity of the 
well is saturated with calcite this drop in aquifer pressure may result 
in the release of C0 2 from the ground water which will then become 
supersaturated. Waters which are supersaturated tend to be unstable 
resulting in the precipitation of carbonates. Some level of super- 
saturation is desirable as the precipitation of carbonates coats the 
distribution system with a protective film reducing corrosion. Excessive 
supersaturation however can result in the encrustation of well screens 
and the precipitation of carbonates in the pore spaces of surrounding 
aquifers reducing well yields. 



Water Quality Characteristics 

The chemical constituents of ground and surface waters are largely 
derived from the solution of minerals in the overburden and bedrock 
through which the waters pass. Both the overburden which is largely 
derived from the local limestone bedrock and the bedrock itself, act as 
sources of carbonate minerals. The ground and surface waters of the 
Wilton Creek drainage basin possess certain chemical characteristics 
that are typical of waters in a carbonate terrain. As rain water 
percolates through the soil horizon, it dissolves large amounts of 
carbon dioxide forming carbonic acid. Carbonic acid is a common solvent 
of carbonates. Analyses of waters in a carbonate terrain are characterized 
by high concentrations of calcium, magnesium and bicarbonate ions resulting 
from the solution of calcite and dolomite. Minor chemical constituents 
in the waters result from the solution of other minerals in the bedrock 
such as gypsum, clay minerals and salts. 

The three major anions, bicarbonate (HCO ) , sulphate (SO ) and 
chloride (Cl~) and the four cations, calcium TCa ++ ) , magnesium (Mg ++ ) , 
sodium (Na + ) and potassium (K + ) are examined in this study. The chemical 
analyses of the water samples are listed in Table 13. The predominant 
cations and anions of the ground-water samples are presented diagram- 
matically in a Piper trilinear plot (Figure 13) . Each vertex represents 
the 100 percent concentration level of a selected ion or group of ions. 
The distribution of the major anions and cations is dicussed in the next 
section. 

The nitrate, phosphorus and iron concentrations in both the surface 
and the ground waters vary considerably over the basin. High nitrate 
and phosphorus concentrations have often been associated with well 
contamination from surface sources such as barnyard wastes or septic 
systems. These generally occur in a very localized manner. High iron 
concentrations in the MOE observation wells are attributed to the 
leaching of iron from the well casing into the surrounding ground 
water. 

Correlation of the Hydrochemistry with the Direction of Ground Water Flow 

Ground water, as it moves through the basin following infiltration, 
is subject to a series of reactive chemical processes which change the 
chemical composition of the water. These processes are discussed in 
considerable detail in Chebotarev (1955), Schoeller (1959) and Hem 
(1970). When viewed on a regional scale, certain sequential patterns 
can be established from hydrochemical data. One such pattern, the 
evolutionary sequence in chemical composition or "metamorphism" as 
described by Chebotarev (1955) for ground waters within a regional flow 
system, is as follows: 

HC0~-*HC0" + C1~-*C1~ + HC0~-*C1~ + S0~ or S0~ + C1~-+C1~ 
3 3 3 4 4 

This is normally accompanied by an increase in the total minerali- 
zation of the water. The application of Chebotarev" s sequence to delineate 
the pattern of ground-water movement in a specific area should be used 
in unison with a previous understanding of the flow system(s). For 
example, individual water samples tend to reflect localized conditions 
which in a heterogeneous hydrogeological environment may be anomalous 
with respect to the sample's position in the flow system. The statistical 
treatment of several such samples may be helpful in reflecting regional 
and possibly basin-wide trends in hydrochemistry. 
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PERCENTAGE OF CONSTITUENTS 

Nomenclature for hydrochemical facies 

(After BACK. 1961 1 




Figure 13 Piper trilinear plot o* chemical analyse* 
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In order to identify and assess hydrochemical trends, the water 
samples are grouped by: (a) hydrochemical similarities; (b) the sampling 
area. The classification scheme proposed by Back (1961) and presented 
in Table 14, is used in this study to establish hydrochemical similarities, 
Water samples are classified into hydrochemical facies based on the 
percentage of predominant anions and cations in the water. This grouping 
is illustrated in the parallelograms which accompany Figure 13. 

TABLE 14. CLASSIFICATION OF HYDROCHEMICAL FACIES (AFTER BACK, 1961) 



Constituent (present) 



Ca+M Na+K HCO ,+CO, S0.+C1 
g 3 3 4 



Cation facies: 

Calcium-magnesium 90-100 0<10 

Calcium-sodium 50-90 10<50 

Sodium-calcium 10-50 50<90 

Sodium 0-10 90-100 

Anion facies: 

Bicarbonate 

Bicarbonate-chloride- 
sulphate 

Chloride-sulphate- 
bicarbonate 

Chloride- sulphate 



90-100 


0<10 


50-90 


10<50 


10-50 


50<90 


0-10 


90-100 



The majority of the ground-water samples fall into the calcium- 
magnesium and calcium- sodium cation facies and the bicarbonate- chloride- 
sulphate anion facies (Figure 13) . These samples were collected from 
wells situated in different areas in the basin and completed to different 
depths. A total of 10 samples fall into the sodium and sodium-calcium 
cation facies, while 9 samples fall into the chloride-sulphate and 
chloride-sulphate-bicarbonate anion facies. These samples were generally 
collected from the deeper wells in the basin. This is a common phenomenon 
in this part of Ontario. Many of the deeper wells in the Belleville- 
Kingston-Cornwall area have been abandoned because high sodium chloride 
concentrations in the well waters have made the waters unpotable. 
Sulphate-type waters with sulphate concentrations in excess of bicarbonate 
concentrations occur in the deeper wells (80 feet+) in the Violet area. 
High sulphate concentrations occur in other parts of the basin; however, 
these waters are generally of the bicarbonate or chloride types. 

In order to establish the distribution of the chemical constituents 
in the ground water in the basin and to illustrate changes in the proportion 
of ions from one area to the next, the sample analyses were grouped by 
the area in which they were collected. These groupings include samples 
collected near the head of the basin - the Harrowsmith - Hartington 
area, the center of the basin - the Violet area, and near the mouth of 
the basin - the Chambers area. As previously indicated ground-water 
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flow is in a downslope direction from upland areas to the valleys and 
from the basin head and divide to the mouth. The results of the stati- 
stical treatment of the analyses listed in these groupings are presented 
in Table 15 along with ion ratios and the calculated levels of saturation. 

The Harrowsmith - Hartington area covers the northern quarter of 
the basin, an area characterized by its strong topographic relief, which 
from the basin divide to the creek valley is in the order of 100-125 
feet. The incised creek valley has well defined bedrock terraces which 
are for the most part obscured by sand and gravel deposits. Ground- 
water flow is from areas of high topographic relief downslope directly 
into the creek. Seepages are common along the valley walls. Water-well 
measurements taken in WC-1 (Figure 11) illustrate a vertical ground- 
water flow vector. 

A total of 22 ground and 4 surface-water samples were collected in 
this area (Map 6, Table 13) . Two of the ground-water samples were 
collected from observation well WC-1. This piezometer nest is estab- 
lished in the sand and gravel plug which acts as the upper basin divide. 
The waters from the two piezometers (25 and 63 feet in depth) have the 
lowest mineral content of any samples collected in the basin. The total 
dissolved solids content for the two samples are 370 and 110 mg/1, 
respectively. The rest of the ground-water samples were collected along 
Highway 38, in order to provide a hydrochemical section through the 
area. Figure 14 is a schematic cross-section illustrating the sample 
location, depth to water (based on well-log information) and the chloride 
and total dissolved solids (TDS) concentrations in the samples. No 
discernable trend could be established from this cross-section. 

The samples collected in the Harrowsmith-Hartington area are 
typically of the calcium-bicarbonate type (Figure 13) with mean HCO /SO , 
HC0 3 /C1 and S0 4 /C1 values of 6.2, 4.0 and 0.9, respectively. 

The Violet area is in the south-central portion of the basin. This 
area is a limestone plain with low relief and little overburden cover. 
The creek dissects the centre of the plain along a fault-controlled 
channel (Figure 3). Relief is in the order of 30 to 50 feet from the 
basin divide to the creek. Ground-water flow is in the direction of the 
creek. 

A total of 17 ground-water and 2 surface water samples were collected in 
the Violet area (Map 6, Table 13). Detailed sampling was conducted 
along a traverse perpendicular to the creek. At this location the MOE 
has installed six observation wells and piezometer nests across the 
bedrock valley. Samples were collected from each of these installations 
and from domestic water wells along the section. The piezometric levels 
in the observation wells were measured prior to sampling. The measured 
piezometric levels and the chloride and TDS concentrations in the ground 
water are shown in Figure 15, a schematic cross-section through this 
area. Ground-water movement along this section is towards the stream 
with a downward gradient on the highs and an upward gradient in the 
valley. As is evident from the piezometer nests WC-3 and WC-6a, the 
chloride and TDS concentrations are higher for the deeper piezometers. 
The upward flow gradient and the high chloride and TDS concentrations in 
WC-3 at depth may indicate the presence of two distinct flow systems, a 
shallow local and an underlying regional system. Chloride and TDS 
concentrations are also higher in the shallow piezometers in the vicinity 
of the creek as opposed to the surrounding uplands. This would in part 
demonstrate an increase in the mineralization of the ground water as it 
moves towards the stream and possibly a mixing phenomenon of the shallow 
and deeper waters. 
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Figure 14. Cross-section along Highway 38 illustrating the chloride and total dissolved solids content in the ground water. 
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Figure 1 6. A north-south cross-section through the observation walla near Violet illustrating the piszometric level* and tha chloride and total dissolved 
solids content in tha ground water. 



For purposes of comparison, the water sample analyses for the 
Violet area are separated into two groups on the basis of hydrochemical 
simlarities. Group A samples are classified as being of the calcium- 
bicarbonate type. These samples were collected from wells ranging in 
depth from 13 to 103 feet. The hydrochemistry of these samples is 
similar to that of the samples collected in the Harrowsmith-Hartington. 
The mean HCO3/SO4, HCO3/CI and SO4/CI values are 5.5, 3.4 and 0.65, 
respectively. 

Group B samples are classified as being of the sodium- sulphate- 
chloride type. These samples are typically more mineralized than those 
of either Group A or of the Harrowsmith - Hartington area. Three of the 
samples WC-3a, WC-3b and WC-6c have HCO3/SO4 values (0.2, 0.3 and 0.1, 
respectively) which are well below those from other areas. Samples WC- 
3a, and DW-29 have chloride contents (2600 and 636 mg/1) which are above 
those of the Harrowsmith - Hartington area and the Group A samples. The 
Group B samples range in depth from 85 to 178 feet. The mean HCO3/SO4, 
HCO3/CI and SO4/CI values are 1.2, 0.3 and 0.6, respectively. 

The Chambers area is in the southern reaches of the basin. This 
area is a near flat plain with relief in the order of 25 feet. Overburden 
deposits such as drumlins, have modified the local topography. Ground- 
water flow in this area has a downward and lateral vector near the basin 
divide and an upward vector in the vicinity of the creek and creek 
mouth. This is borne out in observation well WC-8 where there is a 
strong upward vector to the ground-water flow (Figure 12). 

A total of 10 ground and 4 surface-water samples were collected in 
this area (Map 6, Table 13) . The ground-water samples are divided into 
two groups on the basis of their hydrochemical character. Group C 
samples were collected from shallow (19-60 ft.) wells in this area. 
These wells are typically of the calcium-bicarbonate type. The mean 
HC03/S0 3 , HCO3/CI and SO4/CI values for the samples which fall into this 
category are 14.8, 8.3 and 1.2, respectively. 

Group D, represented by two samples, is characterized by waters of 
the sodium - bicarbonate - chloride type. Two samples were collected 
from observation well WC-8, one from the deep (178 feet) and one from 
the shallow (29 feet) piezometers. Both samples have a high pH measure- 
ment and similar Ca (330 and 320 mg/1) and Mg (75 and 59 mg/1) contents. 
The shallow sample has a distinctive calcium-bicarbonate character, 
whereas the deep samples is of the sodium- chloride type. The mean HCO3/ 
S0 4 , HCO3/CI and SO4/CI values for WC-8f and DW-21 are 9.8, 0.9 and 0.1, 
respectively . 

The following conclusions can be drawn from the distribution of the 
chemical constituents in the ground waters of the basin: 

1. The chemical composition of the ground water changes as the 
water moves through the basin. This change is apparent from 
the grouping of the samples by area. The total dissolved 
solids concentration and the proportion of chloride and sulphate 
concentrations with respect to the bicarbonate concentration 
increase in the direction of ground-water movement. 

2. The chemical concentration (as represented by the TDS concentration) 
in the ground waters is greater for deep wells in the Violet 

and Chambers area then shallow wells. 
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4. 



3. The chemical concentration (based on wells in the Violet area) 
is greater in ground-water samples collected from wells in the 
vicinity of the creek then in samples from wells in the adjacent 
upland areas. Ground-water flow is towards the creek. 
The chemical composition of the ground-water samples collected 
in the upper reaches of the basin and most of the shallow 
wells in the other areas have a similar calcium-bicarbonate 
character. 
The surface water samples collected in the basin are of the calcium - 
bicarbonate type. There is a general increase in the total dissolved 
solids, sodium and chloride concentrations in the water from the north 
to south. The other major constituents in the water are constant. 

Suitability of Ground and Surface Waters for Use 

The chemical quality of the ground and surface waters of an area is 
an important consideration in evaluating the suitability of these waters 
for specific uses. These waters should not contain constituents in 
concentrations which may render them harmful to the use to which they 
may be put. The Ontario Ministry of the Environment has established a 
series of water quality criteria for the Province based on the major 
water use categories. These criteria are presented and discussed in the 
MOE booklet entitled "Guidelines and Criteria for Water Quality Management 
in Ontario" (1974) . The published criteria are listed under permissible 
and desirable concentrations. Waters which exceed the permissible 
concentration criteria should not be used for supply purposes. Good 
quality waters are those which meet the desirable concentration criteria. 
The permissible and desirable concentration criteria for some of the 
major chemical constituents for public water supples are presented 
below: 



Chemical 
Constituents 

Sulphate 

Total Dissolved Solids 

Chloride 

Chloride 

Hardness 

Nitrate as N 
Iron (filterable) 
pH range 



Permissible 
Criteria (in mg/1) 



Desirable 
Criteria (in mg/1) 

<50 



250 

500 

250 <25 

Acceptable levels will vary with local 
hydrogeologic conditions and consumer 
acceptance. 

10 virtually absent 

03 virtually absent 



6.0-8.5 units 



least amount of 
interference with 
treatment process. 



Pollution Indicator 
Organisms 



Coliform and other pollution organisms should 
be virtually absent from all supplies. 



The ground and surface waters in the basin are typically hard (30- 
1850 mg/1 for ground water, 150-322 mg/1 for surface water) . Acceptable 
levels of hardness vary with the use to which the water is put. Of the 
total of 49 ground-water analyses presented in this study, 28 (57%) 
exceed the permissible criteria standards for TDS, 6 (12%) for chloride 
and 1 (2%) for sulphate. The surface water samples meet the permissible 
criteria standards. 

The microbiological content of the water supply should be established 
prior to use by submitting a sample of the water to the local Health 
Unit for examination. Nuisance organisms should be virtually absent 
from water supples. Where they occur it will be necessary to disinfect 
the water supply. 

The suitability of the water for irrigation should be established 
through the local agricultural representative. The relative tolerance 
of crop plants to different chemical constituents in the water varies 
considerable and should be established prior to the initiation of irrigation. 
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SUMMARY 

This report presents the findings of a study conducted on the 
geology and ground-water resources of the Wilton Creek drainage basin. 
The Wilton Creek drainage basin in one of the seven basins studied by 
the Ontario Ministry of the Environment, as its contribution to the 
International Hydrological Decade (IHD) Representative Basin Program. 
Studies conducted on this basin consisted of the collection and inter- 
pretation of geological and ground and surface-water data through field 
surveys, the installation and periodic measurement of eight observation- 
well installations and two streamflow recording stations, and surface 
and ground-water sampling. 

The drainage basin lies within the Napanee Plain physiographic 
region, a flat plain of limestone which dips gently to the southwest 
resulting in a series of shallow erosional scarps. Surface relief is 
provided by drainage and glacial erosion. The overburden is thin and 
discontinuous for the most part; however, local thicknesses in the order 
of 90 feet have been recorded in the southern reaches of the basin. 

The climate is influenced by the moderating effects of Lake Ontario. 
The mean annual temperature at Hartington for the interval 1968 to 1973 
is 43.3°F. The mean annual precipitation for that interval is 36.33 
inches at Hartington and 36.44 inches at Violet. 

The bedrock is comprised of Precambrian metamorphic and igneous 
rocks overlain non-conformably by Cambrian and Ordovician sedimentary 
rocks (Map 1). The sedimentary succession consists of conglomerates, 
sandstones, shales and carbonates of the Basal and Simcoe groups. The 
Basal Group includes the Potsdam and Shadow Lake formations. The 
Simcoe Group includes four rock units, the Lindsay, Verulam, Bobcaygeon 
and Gull River formations. The uppermost unit, the Lindsay Formation, 
is not present in the basin area. 

The Pleistocene overburden includes glacial tills, glaciof luvial 
sands and gravels and stratified glaciolacustrine silts and clays (Map 
2). Deposition appears to have occurred along an oscillating ice front. 
The lack of Lake Iroquois beach and near-shore deposits suggest that the 
ice margin may have been floated clear of the land surface by the buoyant 
action of Lake Iroquois waters. A single lacustrine shore feature, a 
gravel bar south of Morven, has been related to a brief transitory post- 
Lake Iroquois stage, Lake Belleville. 

Ground water is the major source of domestic water supply. The 
carbonate bedrock of the Simcoe Group constitutes the principal aquifer. 
Ground-water availability and quantity is related to the number of 
interstitial openings such as fractures in the aquifer and their degree 
of interconnection. The density of openings varies both with location 
and depth in the aquifer. 

The calculated coefficients of transmissibility range from 90-3900 
Igpd/ft. A single permeability value is in the order of 100 Igpd/ft 2 
and the coefficient of storage is estimated to be 4 x 10 -4 . The mean 
coefficient of transmissibility as estimated from specific-capacity 
values is 2,300 Igpd/ft. 



On the basis of productivity- frequency graphs prepared from adjusted 
specific-capacity values and the statistical treatment of these values, 
significant relationships exist between the well productivity and 1) the 
depth of penetration by the well, 2) the position of the well relative 
to topography and 3) the bedrock units encountered. Shallow wells and 
wells finished in valleys are more productive than deep wells and wells 
finished in upland areas. 

The till and lacustrine deposits act as a low permeable mantle, 
creating semi-confined aquifer conditions. The sand and gravel deposits 
constitute locally favourable aquifers; however, their thickness and 
areal extent are too small to yield anything but small domestic water 
supplies. The mean coefficient of transmissibility of the overburden, 
as estimated from specific-capacity data, is 1400 Igpd/ft. The mean 
permeability as estimated from grain-size distribution is in the order 
of 860 Igpd/ft 2 . 

Ground-water flow in the basin is strongly influenced by the local 
topographic relief. Ground water moves from upland areas to valleys and 
in a more general manner from the head of the basin to the mouth. The 
ground-water level fluctuates seasonally with the greatest rise occurring 
during the spring-melt period. Water-level fluctuations as measured in 
observations wells vary from 0.7 to 34.0 feet over the basin, depending 
on well location and depth. There are no noticable long-term changes in 
ground-water levels over the period 1967 to 1973. 

A simplified hydrologic budget was estimated for the water-year 
interval 1967-68 to 1973-74 using the following equation: 



P = R„ + R„ + ET 



where: 



P = precipitation 

Rg = ground water runoff or baseflow 

R s = surface runoff 

ET = evapotranspiration 

Precipitation for the seven-year period, as averaged for Hartington 
and Morven, is 37.04 inches/year. The mean total streamflow runoff for 
the same period, as measured at hydrometric station 3MH-004, averaged 
16.40 inches/year, while the ground-water runoff component (estimated 
from the monthly minimum streamflow) averaged 3.97 inches/year (24.2% of 
total streamflow) . The potential evapotranspiration as estimated by the 
Munson-Index method is 22.68 inches/year. The evapotranspiration loss 
calculated by the water-balance equation is 20.64 inchs/year. The 
difference (2.04 inches/year) between the potential evapotranspiration 
and that calculated by the water-balance equation represents the cumulated 
errors inherent in the estimation of the components of the budget and 
net changes in basin storage and basin underflow, which have not been 
measured directly. 
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A total of 38 ground-water and 15 surface water samples were 
collected at different locations in the basin to establish the chemical 
composition of the waters. The ground-water samples were arbitrarily 
divided into three groups based on location and hydrochemical similarities 
The following trends became apparent by this grouping: 

1. The chemical composition of the ground water changes as the 
waters move through the basin; TDS content increases in a 
downgradient direction. 

2. The TDS concentration is greater in water samples collected 
from deep wells in the Violet and Chambers area than shallow 
wells. 

3. Ground water in the vicinity of the creek has a greater mineral 
content then that in the adjacent upland areas. 

4. Ground water collected in the upper reaches of the basin and 
in most shallow wells, is of similar composition, being of the 
calcium-bicarbonate type. 
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APPENDIX A 

GUIDE TO THE BEDROCK 

AND OVERBURDEN GEOLOGY 

OF THE 

WILTON CREEK DRAINAGE BASIN 

(compiled during 1973) 

Table 16. Stratigraphic units present at site localities A to I. 

Figure 16. Location map for sites referred to in the Guide. 

Figure 17. General diagram of the landfill site near Violet with 
locations of the sections cited in the Guide. 
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TABLE 16. STRATIGRAPHIC UNITS PRESENT AT SITE LOCALITIES A TO I 



Stratigraphic Sequence 
(Liberty, 1971) 



Relative Stratigraphic Positions of Localities 
A through I 



Verulam Formation 



Bobcaygeon Fm 



UppeT 



Lower 



Member D 



Member C 



Member B 



Upper submember 



i i 



r 



Middle submember 
Lower submember 



Shadow Lake Formation 



E F ! 



Potsdam Formation 



Precambnan 
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Figure 16. Location map for sites referred to in the Guide. 
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Guide to the bedrock geology of the Wilton Creek drainage basin. 
Stratigraphic distribution is presented in Table 16 and site locations 
in Figure 16. 

A. PRECAMBRIAN 

The Precambrian is exposed in a sand and gravel pit in the Wilton 
Creek valley, northeast of Harrowsmith. 

The Precambrian at this locality is composed of granite gneiss in 
the form of an inlier with an estimated relief of about 75 feet 
(determined from adjacent water-well logs) . 

Overlying the Precambrian unconformably , are the Shadow Lake and 
Gull River formations (Text Photo la) . 

Sequence Exposed 



Gull River Member A 

Shadow Lake 
Precambrian 



upper submember 
middle submember 



The upper submember of Member A may be observed higher up on the 
valley wall east of this location, at approximately 550 feet. The 
lower 20 feet of the Gull River Formation is not exposed at this 
location. 

The absence of the Potsdam Formation is explained by the onlapping 
of the Paleozoic sequence over the Precambrian highs. The Potsdam 
Formation although occurring near the base of the inlier, thins out 
higher on the flanks (see Text Figure 2) . 

B. POTSDAM FORMATION 

Although the Potsdam Formation does not outcrop in the basin, it is 

inferred to be present from water-well logs. 

The Potsdam Formation examined to the northeast of the basin on the 

shore of Knowlton Lake, consists of a fine to medium - grained red 

quartz sandstone. The quartz grains are well-rounded. Hematite 

and calcite make up the cementing matrix. 

Numerous exposures of the Potsdam Formation occur north of Kingston 

(Text photos 2 and 3) . 

C- GULL RIVER FORMATION (MIDDLE AND UPPER SUB-MEMBERS OF MEMBER A) 



These submembers are exposed at a roadcut about 1/2 mile north of 

Harrowsmith on Highway 38. 

The contact between the submembers is represented by a 2-foot thick 

layer of green-black shale. 

The elevation for the contact is approximately 500 feet. 

There appears to be little lithological variation between the two 

submembers, although the middle submember is thicker. 

Beds of the middle submember, about 10-15 feet below the contact, 



are fossilferous and can be used as marker beds locally. 






D. GULL RIVER FORMATION (CONTACT BETWEEN MEMBERS A AND B) 

The contact between the upper submember of Member A with Member B 
is exposed about 3/4 mile north of Murvale on Highway 38. The 
contact is marked by colour changes in the limestone. 
The lithology is repetitive with the exception of a decreasing 
calcite vug content above the contact. 

E. GULL RIVER FORMATION (MEMBER B) 

The member is partially exposed 3/4 mile south of Harrowsmith on 

Highway 38 . 

It consists of alternate beds of massive limestone and thin beds of 

friable and argillaceous limestone which is fossiliferous and 

contains numerous calcite vugs. 

The exposure appears to be the upper submember. 

F. GULL RIVER FORMATION (MEMBER B) 



Member B is partially exposed west of Harrowsmith in a series of 
outcropping erosional scarps made up of hard beds, while the soft 
beds are not exposed. 

GULL RIVER FORMATION (CONTACT BETWEEN MEMBERS B AND C) 

This bedrock exposure occurs as a terraced surface in a sand and 

gravel pit southeast of Wilton. 

The contact between Members B and C is exposed on the bedrock 

terrace at approximately the 450-feet contour. 

Member C is more "arkosic" and is highly fossiliferous. 

Stromatolites are exposed near the abandoned office building of the 

sand and gravel pit. 

GULL RIVER FORMATION - BOBCAYGEON FORMATION CONTACT 

This contact is visible on Highway 401 in a road-cut east of the 
Highway 401 and Highway 133 intersection. 

The contact is distinctive between the more massive Gull River 
(Member D) and the friable argillaceous Bobcaygeon limestone (Text 
Photo 11) . 

GULL RIVER - BOBCAYGEON - VERULAM FORMATION 

These formations are all exposed in the vicinity of Chambers. 

The Gull River Formation, Member D, is visible in the streambed as 

is the Bobcaygeon, while the Verulam Formation is exposed further 

south in two abandoned quarries. 

Text Photo 12 shows the Verulam Formation. 
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GUIDE TO OVERBURDEN GEOLOGY OF THE WILTON CREEK DRAINAGE BASIN 

1. Sand, Gravel and Till "Plug" ■ 

At this location, a sand, gravel and till plug fills the Wilton 
Creek Valley and acts as a drainage divide between the Creek and 
Knowlton Lake. 

The plug appears to act as the source area for Wilton Creek. 
The following section is exposed in a small pit on the Knowlton 
Lake side of the plug: 

Thickne_ss_ 

2-10 Feet: Sandy till with few stones, very similar to the under- 
lying glaciofluvial sands (Kame deposit ?) . 

20 + Feet: Medium-grained sand, good sorting. 

2 - Sand and Gravel Pit (1-1/2 miles east of Hartington) 

At this locality sand and gravel fills the base of the valley and 

occurs along its terraced bedrock sides. 

This sand and gravel appears to be a continuation of the plug down 

the creek valley. 

The exposed section is presented below: 

Thickness 

2 Feet: Sandy till with few stones, local thickness to 5 feet. 
15 Feet: Medium sand, well sorted. 

5 Feet: Coarse sand and gravel, layer of well-rounded boulders 
which is assumed to be ice contact in origin. 

3'« Kame Terrace Sand and Gravel and a Limestone Quarry 
(1 mile east of Hartington) 

The limestone quarry is in the middle submember of Member A of the 

Gull River Formation. 

The sand and gravel along the side of the valley is deposited on 

bedrock terraces. 

The exposed section is presented below: 

Thickness 

2 Feet: Sandy till with few boulders. 

3-5 Feet: Medium sand, stratified and well sorted. 

10 Feet: Coarse sand and gravel, poor to fair sorting. 

There is water seepage from the base of the coarse sand. 



4. Sand and Gravel Pit (northeast of Harrowsmith) 

Kame terrace deposit which has been largely exploited. 

No undisturbed sections present. 

Sand and gravel can be found on a series of terraces up to an 

elevation of 575 feet. 

5. Active Sand and Gravel Pit (near Harrowsmith) 

Sand and gravel deposited along the side of the creek valley in the 

form of a kame terrace deposit. 

The bedrock is exposed in several areas. There are numerous fractures 

and glacial striations. 

Fracture and striae trends are presented below: 

Fracture Trends : 

017°. 357°, 347°, 334°, 012°, 020°, 014°, 012° 

Striations: 

234°, 237°, 242°, 240° 

The exposed section at this locality is as follows: 

Thickness 

2-5 Feet: Sandy till with few stones, strongly weathered. 
10-20 Feet: Medium sand, cross bedded with forsets dipping to the 
southeast, sand rests on bedrock. 

Sand thins out towards the east and the south, local pockets of 
coarser sand and gravel. 
Seepage at base of the sand. 

6. Abandoned Sand and Gravel Pit (two mii^s southwest of Harrowsmith 
on west side of Highway 4-A) 

Kame terrace deposit along the side of the valley. 
- Exposed section is presented below: 

Thickness 

25 Feet: Sandy till with few stones. 

<50 Feet: Sand and gravel has been removed; thickness is assumed. 

7. Active Sand Pit (north of Wilton on Highway 4-A) 

This appears to be a continuation of the sand terraces found both 

to the southwest and northeast of this deposit. 

Deposit is exposed on both sides of the road, (with sand resting on 

the bedrock terraces) . 

The section on the south side of the road is presented below: 
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Thickness 

3 Feet: Sandy till with few boulders, soil profile developed. 

1 Foot: Stoney till, dense. 

10 Feet +: Fine to medium sand, stratified and well sorted. 

Ground-water seepage is evident at the base of the sand. 
Fracture and striae trends are presented below: 

Fracture Trends: 

347°, 357°, 358°, 357°, 358° 

Striations: 

222°, 227° 

The section on the north side of the road is presented below: 
Thickness 

1 Foot: Till (upper surface has been stripped so thickness is not 

measureable) . 
15-20 Feet: Medium-grained sands, foreset beds dip to the southwest. 

8. Sand and Gravel Pit (south of Wilton on east side of creek) 

Most of the sand and gravel has been removed from this terrace 

deposit. 

No evidence of till. 

Water table is near the base of the sand. 

The section is presented below: 

Thickness 

4 Feet: Silty clay layer, strongly weathered. 

1/2 Foot: Massive dense clay. 

10 Feet +: Sand and silt, medium-grained sand is well stratified, 

coarser-grained sand underlies the medium-grained sand 

locally. 

9 - Sand and Gravel Pit (southwest of (8) about 1/2 mile) 

The exposed section is presented below: 
Thickness : 

4-5 Feet: Sandy-silty till, well developed soil horizon. 
25-35 Feet: Medium, cross-bedded sand, foreset beds dip to southwest. 
10-15 Feet: Coarse sand and gravel of which the lower 5-8 Feet are 
saturated. 



The pit must be dewatered continuously during the spring and early 
summer (Photo 24) . 

10. Large Drumlin (southeast of Thorpe) 

Drumlin has a stratified sand and gravel core covered by a 2 to 5- 
foot till layer. 

Sand and gravel pit developed on the northeast end of the drumlin 
with about 20 feet of unsorted sand and gravel exposed (Photo 16) . 

11. Sand and Gravel Pit (northeast of the intersection of h ighways 133 
and 401) . 

Well developed bedrock terraces upon which interbedded clay, silt, 

sand and gravel have been deposited. 

Sand is stratified, foreset beds dip to southwest. 

Measured sections are presented below: 

Southeast End Of Pit ; 

Thickness 

4 Feet: Sandy, silty till, few boulders, weathered profile. 

gradational contact. 

4 Feet: Medium, well-sorted sand. 

Sharp contact. 

1 Inch: Dense clay. 

9 Inches: Interbedded sand and reddish brown silty-sand till. 

1 Foot: Medium-grained sand. 

2 Feet: Silt and clay. 
2-1/2 Feet: Medium-grained sand, becomes increasingly siltier towards 

the base. 



Scalloped contact. 

1 Foot: Silty till with clay layer (1 inch). 

4 Feet: Medium sand and gravel near base. 

5 Inches: Dense red-brown silty till (text Photo 14) clay, sand, 

till, red clay, till, sand. 
9 Inches: Clay and sand. 
4 Feet: Medium sand. 
1 Inch: Gravel. 
1 Foot: Medium and coarse sand. 

Sand lies on a limestone terrace. 

Central portion of the pit may contain another 10+ feet of sand on 

a lower terrace; however, it was flooded during field work. 

Northeast End of the Pit : 

Thickness 

3 Feet: Silt and sand (possible a strongly weathered till) . 
3 Feet: Medium sand. 

8 Feet: Stoney clay till interbedded with sands above a bedrock 
terrace. 
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Figure 17. General diagram of the landfill site near Violet with locations of the sections 
cited in the Guide. 
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There is considerable seepage at the base of the sand. 
Fracture trends and striae measurements are given below: 

Fracture Trends: 

030°, 040°, 300°, 042°, 300°, 042° 

Striations: 

o o o 
236 , 234 , 238 

12. Sanitary Landfill Site (near Violet) 

Originally a sand and gravel pit, presently being backfilled with 
solid waste. 

There is ground-water seepage from the base of the sand near the 
bedrock contact. 



Fracture Trends: 

o 

, 32? 
.o 



320°, 328 , 022°, 292°, 299°, 289°, 004°, 354°, 339°, 344°, 277°, 277°, 
287 C 



Striations: 



o o 
236 , 234 



Various measured sections A - F are presented as follows: 
Section A. 



Thickness 

3 Feet: 

2 Feet: 

3 Inches: 
6 Inches: 
8 Inches: 
18 Feet: 

Section B. 

Thickness 

3 Feet: 
6 Inches: 
3 Inches: 
6 Inches: 

1 Foot: 

2 Feet: 
18 Feet: 



Weathered silt - clay layer with soil profile. 

Cross-bedded silt and clay - reworked clay fragments at 

the base of this layer. 

Till - dense red silty-clay till, abundant stones. 

Medium sand - sand grades into till. 

Till - stoney, sand. 

Medium sand, good sorting, base is not exposed. 



Weathered silty clay, soil profile. 

Clay till, moderate stone content. 

Stratified sand and silt layer. 

Compact silty till. 

Coarse sand. 

Gravel foreset beds dip the southwest 

Medium and fine sand. 
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Section C. 



Up to 20 feet of coarse sand. 

Location is on a lower bedrock terrace, the base of which is flooded. 
Sand and gravel can be traced along bedrock terrace southwest to 
County Road. 



Section D. 
Thickness 

5 Feet: 
2 Feet: 

6 Inches: 
9 Inches: 

1 Inch: 

7 Inches: 

2 Inches: 
20 Feet: 

Section E 
Thickness 

4 Feet: 

6 Inches: 

5 Feet +: 



Silty clay, soil profile. 

Brown clay, evenly-bedded, graded bedding. 

Silty clay. 

Stoney sandy till. 

Dense, plastic clay, red-brown colour. 

Stoney sandy till. 

Gravel layer. 

Medium, cross-bedded sand - undulating (scalloped) surface 

suggest overriding. 



Silty clay, evenly-bedded, soil profile. 
Stoney till, dense, some clay balls. 
Clay-silt varves, good bedding. 



Base is buried. It is assumed that till underlies the clay, as 
inferred from adjacent area. 

Varves slightly west of Section E appear to represent a local ponding in 
the immediate area. 

13 • "Beach Remnant" (north side of drumlin near the borders of the 

townships of North Fredericksburg, South Fredericksburg and Ernestown) 

This site, described by Mirynech (1967), is a wave-modified Quinte 
Bay drumlin. 

This "beach remnant" suggests a brief water stand, the Belleville 

Stage (Mirynech, 1962) . 

A wave-deposited gravel bar present on the northeast flank of the 

drumlin extends for about 1 mile. 

The deposit has been exploited for the gravel. 

14 • Qunite Bay Drumlin (on Highway 8 near the lower edge of the basin) . 

Typical long, low-lying drumlin about 30-40 feet high. 
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APPENDIX B 

PUMPING TEST DATA 

Figure 18. Time-drawdown curve of an observation well near WC-7. 

Figure 19. Time-drawdown data for pumped well WC-7. 

Figure 20. Time-drawdown data for pumped well WC-4. 

Figure 21. Time-drawdown data for pumped well at service centre 
E-5. 
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Field data 



r Nonequilibnum type curve 

T _ 114 6 QW (u ' / B) 
S 
where T -coefficient of transmissibility (imperial gallons day /ft) 

Q=- discharge of the well (3 imperial gallons /minute) 

s -drawdown (0 68 feet) 

W (u, '/ B) - well function for leaky artesian aquifers 

(2 8 dimensionless). obtained from graphical superposition 

T 114 6(3)(2 8) 
68 

■ 141 5 6 imperial gallons per day per fool 
1 87 r 
where S - coefficient of storage (dimensionless) 

t -time since start of pumping (12 2 minutes) 

r -distance from the discharging well to the 
point of observation (40 feet) 

u -well function ('/(j- 10 0) 



1415 6(10"') (12 2) 1 



1 87 (1600) 
■4 Ox 10 -4 



1440 



(Convert minutes to days) 



(Hantush and Jacob, 1955) 



5 6 7 8910 



Time in minutes 



Figure 1 8. Time-drawdown curve of an observation well near WC-7 
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__ Fitted line 

T= 264Q 

AS 
T= coefficient of ttansmissibility 
Q = pumping late 
AS= change in drawdown 

Q= 3 Imperial gallons per minute 
AS= 8.6 feet 
T= 264(3) 
8.6 
"-92 Imperial gallons per day per foot 

(Cooper and Jacob, 1946) 
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Figure 19. Time-drawdown data for pumped well WC-7. 
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T = coefficient of transmissibillty 
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— j AS- change in drawdown 

Q 3 Imperial gallons per minute 
AS= 9.6 feet 
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Figure 20. Time-drawdown data for pumped well WC-4. 
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T = 264Q 
AS 

T= coefficient of transmissiblllty 

O = pumping rate 
AS= change in drawdown 

Q = 100 Imperial gallons per minute 
AS = 8.7 feet 

T= 264(100) 

6.7 
^3900 Imperial gallons per day per foot 

(Cooper and Jacob. 1946) 
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Figure 21. Time-drawdown data for pumped well at service centre E-5. 
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APPENDIX C 

WATER QUALITY ANALYSES 

Table 13. Results of chemical analyses of water samples. 

Table 15. Grouping of water analyses by their areal distribution in 
the basin. 
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10 

11 

12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
35 



2 tt) 
O H 
M J3 



Birdsell 
Hartington P.S. 
Leonards Garage 
Hartington Farm 
Supply 
R. Goslin 
Hanna 
Bennett 
Harrowsmith 
Restaurant 
Harrowsmith 
Post Office 
Robinson 
Simmon's 
(Lamberts' Gar.) 

C. Hanson 
Wistard 

D. Gowdy 

F. Kent 
Harden 
Dartt 
Steele 

A. Alkenbrack 
L. Wannamaker 
J. Chambers 
P. Wycott 
L. Lasher 

G. Lasher 
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78 
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53 
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60 

70 
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48 
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19 
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65 
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8.0 
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7.3 
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7.9 
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Constituents in Milligrams/litre (mg/1) 



246 
360 
453 
235 

255 
275 
364 
434 

360 

448 
372 
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78 
58 
47 
73 
30 
20 

191 

20 

7 

463 
19 
71 

356 



150 

55 

21 
210 

57 
29 
51 
59 

100 

60 
41 

44 
170 
50 
110 
38 
39 
58 
49 
50 
60 
38 
51 
74 



0.1 
0.1 
0.08 
0.28 

0.1 
1.5 
1.5 
0.1 

5.6 

1.7 
0.1 

0.3 

0.4 

0.5 

2.1 

3.6 

1.1 

0.1 

3.0 

0.1 

9.5 

0.02 

0.1 

0.01 



0.012 
0.001 
0.010 
0.072 

0.002 
0.001 
0.001 
0.028 



800 
660 
830 
850 

460 
470 
700 
740 



0.026 820 



0.050 
0.022 

0.001 
0.001 
0.001 
0.001 
0.004 
0.001 
0.008 
0.037 
0.022 
0.066 
0.002 
0.001 



680 
730 

440 
650 
460 
550 
400 
390 
800 
460 
340 
2130 
340 
450 
960 



1100 

820 

1080 

1000 

600 

640 

910 

1020 

1100 



1300 
1050 
1370 
1340 

750 

780 

1180 

1220 

1300 



950 1140 
850 1040 



550 
780 
620 
720 
520 
540 
890 
520 
490 
3000 
550 
620 



Wa : 



710 
980 
780 
910 
640 
630 

1410 
700 
550 

3350 
530 
840 

1750 

er So 



30* 
436 
532 

138 

246 
328 

252 
472 

484 

496 
480 

282 
360 
336 
376 
328 
328 
102 
348 
266 
576 
246 
270 



DW 2 4 

DW 2 5 

DW 26 

DW 27 

DW 2 8 

DW 29 

DW 30 

DW 3 1 

DW 32 

DW 3 3 

DW 34 



. § 

& IT] 



? Shallow 



14 
15 
18 
50 
60 
90 
60 
70 
75 
70 



July 10/6? 



7 . 8 

7. 4 
7 . 4 
7.6 
7 . 8 

8. 5 
7 .9 
7 . 8 
7 . 7 



326 
225 
269 
385 
386 
513 
273 
224 
213 
281 
459 



Constituents in Milligrams/litre (mg/1) 



O "3 
3° 



166 

80 

109 

142 

143 

162 

90 

80 

77 

90 

150 



46 
14 
16 
10 
10 
45 
32 
13 
24 
27 
25 



28 

13 

17 

37 

43 

452 

37 

11 

7 

5 

53 



< 

e-> x. 
o 



4. 2 

2 .0 
1. 7 
5 . 2 



< <v 

O 



0. 15 
0. 35 
0. 10 
3. 50 



3 


6 


0. 
1 . 2 



2. 

0i 

2i 
li 



398 
274 
328 
470 
471 
626 
333 
27 3 
260 
343 
560 



53 

2 5 

3 
63 
63 

636 
64 
23 
32 
14 
79 



exh 
13 
45 
28 
30 

118 
46 
13 
46 
33 
54 



►J 

< 

E- Z 



exh . 
330 
460 
570 
580 
&670 
550 
330 
400 
415 
730 



w o 



Q U 
55 O 

O -H 
U £ 



e xh . 

115 

700 

896 

922 

[2990 
747 
495 
555 
606 

1110 



*-£> 









St. 



3 
4 

I •: 

5 

6 

7 

1 

12 

e 

9 

10 
11 





s 




r-4 


(--. 






>Tl 




u 




•H 


■' > 




. J 



















«1 




M 








f 



Harrowsmith 



Wilton 



Violet 



Fredricksburgh 



Nov. 
Dec. 
Apr. 
July 
Nov. 
Dec. 
Apr. 
Oct. 
July 
Nov. 
Dec. 
Apr. 
July 



29/73 

11/68 

28,29/69 

10/69 

29/69 

11/68 

28,29/^8 

17/73 

10/69 

29/68 

11/68 

28.29/p9 

10/69 






8.1 
8.0 
8.2 
8.8 
8.1 
8.0 
8.3 
8.5 
8.4 
8.1 
8.2 
8.8 
8.6 



240 



183 
282 
226 
247 
154 
242 
210 
210 
225 
120 
256 
216 
222 



Constituents in Milligrams/litre (mg/1) 



65 
►6 
76 

86 
62 

80 
74 

MM 

SO 

52 
M 
76 

'■-' 



13 

JO 
15 
19 

5 
L7 
10 
16 
13 

5 
14 
12 
L4 



3 

7 

5 

7 

2 

6 

4 

20 

12 

6 

9 

10 

14 



2.1 
2.1 
1.7 
1.4 
2.0 
1.8 
1.7 
4.6 
1.9 
3.1 
2.3 
1.9 
3.0 



0.39 
0.10 
0.15 
0.70 
0.59 
0.30 
0.20 
0.30 
0.20 
1.57 
0.35 
0.25 
0.25 



223 
344 
276 
301 
188 
295 
256 
293 
256 



263 

271 



13 
10 
exh. 

3 

8 

8 
52 
24 
li 
16 
15 
26 



17 

24 
- 
L4 
14 
20 

17 

41 

13 
21 

;: r 5 
20 
20 



0.6 
1.5 
1.6 

0.35 
0.78 
0.04 

0.1 

0.82 
0.72 
0.02 



245 
385 
265 
350 
208 
310 
235 

340 
234 
345 
245 
340 



216 

249 

176 
272 
226 
284 

150 
292 
239 



597 
453 
515 
324 
525 
433 
480 
495 
309 
582 
453 
515 



TABLE 15 



HARROWSMITH-HARTINGTON AREA 



No. Title 



WC1A 
WC1B 
DW01 
DW02 
DW03 
DW04 
DW05 

8 DW06 

9 DW07 

10 DW08 

11 DW09 
DW10 
DW11 

14 DW12 

15 DW13 

16 DW14 

17 DW15 

18 DW16 

19 DW17 

20 DW18 

21 DW33 

22 DW34 



12 
13 



Ca 

77.0 

17.0 

9.0 

125.0 

155.0 

37.0 

66.0 

98.0 

52.0 

123.0 

131.0 

131.0 

149.0 

83.0 

86.0 

88.0 

117.0 

106.0 

92.0 

27.0 

90.0 

150.0 



Mg 

2.0 
1.0 
2.0 
30.0 
35.0 
11.0 
20.0 
20.0 
30.0 
40.0 
38.0 
41.0 
26.0 
18.0 
35.0 
28.0 
20.0 
16.0 
24.0 
8.0 
27.0 
25.0 



Na 

12.0 
5.0 

282.0 
56.0 
84.0 

230.0 
69.0 
37.0 

164.0 
91.0 
91.0 
64.0 
40.0 
39.0 
78.0 
32.0 
39.0 
8.0 
8.0 

276.0 

5.0 

53.0 



Mean 91.3 22.5 80.1 

VIOLET AREA GROUP B SAMPLES 



No. Title 



Ca 



Mg 



DW29 
WC3A 
WC3B 
WC6C 



162.0 45.0 
520.0 132.0 



15.0 
80.0 



8.0 
4.0 



194.3 47.3 



Na 

452.0 

768.0 

99.0 

120.0 

359.8 



6.5 
2.2 
3.2 
4.1 
2.8 
1.3 
4.5 
1.9 

13.0 
7.5 

11.0 
8.0 
2.0 
3.1 

10.0 
4.2 
5.5 
1.7 
2.0 
5.2 
2.4 

30.0 

6.0 



8.6 
15.0 

3.3 
13.0 



CL 

18.0 

5.0 

158.0 

104.0 

175.0 

151.0 

57.0 

70.0 

140.0 

114.0 

145.0 

80.0 

104.0 

78.0 

58.0 

47.0 

73.0 

30.0 

20.0 

191.0 

14.0 

79.0 

86.6 



CL 

636.0 

2600.0 

119.0 

134.0 



10.0 872.3 



ALK 

251.0 
39.0 
246.0 
360.0 
455.0 
235.0 
255.0 
275.0 
364.0 
434.0 
360. C 
448.0 
372.0 
219.0 
271.0 
289.0 
273. C 
255. C 
279.0 
385.0 
281.0 
459. } 

309.3 



AJJC 

513.0 
6.0 

32.0 

31.0 

14E .5 



so 4 

25.0 
35.0 

150.0 
55.0 
21.0 

210.0 
57.0 
29.0 
51.0 
59.0 

100.0 
60.0 
41.0 
44.0 

170.0 
50.0 

110.0 
38.0 
39.0 
58.0 
33.0 
54.0 



S °4 

118.0 

35.0 

110.0 

300.0 

140.7 



NO 

0. 1 
0.1 
0.1 
0.1 
0. 1 
0.3 
0.1 
1.5 
1.5 
0. 1 
5.6 
1.7 
0.1 
0.3 
0.4 
0.5 
2.1 
3.6 
1.1 
0.1 
0.0 
0.0 



67.6 0.9 



Ca/ 

"g 

2.18 

2.39 

1.14 

12.13 

4.4 



Ca/ 

M 9 

23.35 
10.31 
2.73 
2.53 
2.69 
2.04 
2.00 
2.97 
1.05 
1.87 
2.09 
1.94 
3.48 
2.80 
1.49 
1.91 
3.55 
4.02 
2.33 
2.05 
2.02 
3.64 

3.77 



HCO3/ 

S °4 

9.64 
1.07 
1.57 
6.28 
20.79 
1.07 
4.29 
9.10 
6.85 
7.06 
3.46 
7.17 
8.71 
4.78 
1.53 
5.55 
2.38 
6.44 
6.87 
6.37 
8.17 
8.16 

6.24 



HC0 3 / 

S °4 

4.17 
0.16 
0.28 
0.10 

1.20 



HCC3/ 
CL 

0.57 
0.00 
0.19 
0.16 

0.25 



HCO3/ 
CL 

9.88 
5.53 
1.10 
2.45 
1.84 
1.10 
3.17 
2.78 
1.84 
2.70 
1.76 
3.97 
2.53 
1.99 
3.31 
4.36 
2.65 
6.02 
9.88 
1.43 
14.22 
4.12 

4.03 



S0 4 / 
CL 

0.14 
0.01 
0.68 
1.65 

0.64 



SO4/ 
CL 

1.03 
5.17 
0.70 
0.39 
0.09 
1.03 
0.74 
0.31 
0.27 
0.38 
0.51 
0.55 
0.29 
0.42 
2.16 
0.79 
1.11 
0.93 
1.44 
0.22 
1.74 
0.50 

0.94 



ERROR % SATC 



-12.62 

•15.80 

1.80 

-0.15 

-0.46 

-2.14 

0.99 

0.08 

0.55 

1.70 

0.50 

1.25 

0.61 

-0.77 

0.93 

0.10 

-3.29 

0.05 

-0.47 

-0.53 

2.15 

0.36 



ERROR % SATC 



1.64 
-2.48 
-4.10 
-3.80 



0.838 

-1.391 

0.120 

0.630 



2.218 
0.964 
0.426 
0.592 
0.466 
0.247 
0.507 
0.513 
0.323 
0.448 
0.486 
0.391 
0.484 
0.548 
0.491 
0.382 
0.452 
0.424 
0.595 
0.255 
0.688 
0.953 

0.584 



SATD 

1.364 

-3.121 

0.171 

0.161 



SATD 

2.900 
0.859 
0.397 
0.800 
0.530 
0.184 
0.721 
0.566 
0.641 
0.649 
0.670 
0.518 
0.448 
0.659 
0.814 
0.497 
0.363 
0.256 
0.837 
0.214 
1.082 
1.366 



SATG 

-1.477 
-1.768 
-2.162 
-1.159 



SATG 

-2.293 
-2.481 
-2.378 
-1.734 
-2.109 
-1.647 
-1.896 
-2.040 
-2.109 
-1.745 
-1.497 
-1.712 
-1.794 
-1.907 
-1.400 
-1.861 
-1.432 
-1.875 
-1.930 
-2.310 
-2.012 
-1.690 



0.726 -1.902 



SATM 

-0.048 
-2.209 
-0.427 
-0.948 



SATM 

0.204 
-0.584 
-0.508 
-0.271 
-0.414 
-0.541 
-0.264 
-0.425 
-0.160 
-0.278 
-0.295 
-0.351 
-0.515 
-0.368 
-0.155 
-0.363 
-0.567 
-0.646 
-0.237 
-0.519 
-0 . 084 
-0.066 

-0.337 



0.052 -0.362 -1.609 -0.899 






vx> 

ON 



VIOLET AREA GROUP A SAMPLES 



No. Title 



WC2 

WC3C 

WC3D 

WC4 

WC5 

WC6A 

WC6B 

WC7 

DW22 

10 DW2 3 

11 DW30 

12 DW31 

13 DW32 

Mean 



Ca 

85.0 

94.0 

72.0 

200.0 

114.0 

23.0 

15.0 

114.0 

72.0 

81.0 

90.0 

80.0 

77.0 

85.9 



Mg 

17.0 
11.0 
26.0 
24.0 
22.0 
15.0 
7.0 
21.0 
16.0 
16.0 
32.0 
13.0 
24.0 

18.7 



Na 

17.0 
28.0 
14.0 
58.0 
71.0 
93.0 
131.0 
35.0 
15.0 
55.0 
37.0 
11.0 
7.0 

44.0 



K 

2.3 
1.9 
3.9 
4.8 
6.5 
7.6 
9.8 
3.8 
2.9 
4.0 
4.0 
1.9 
3.4 

4.3 



CL 

37.0 

61.0 

35.0 

255.0 

133.0 

71.0 

85.0 

126.0 

19.0 

71.0 

64.0 

23.0 

32.0 

77.8 



ALK 

227.0 
206.0 
217.0 
268.0 
290.0 
190.0 
195.0 
247.0 
222.0 
247.0 
273.0 
224.0 
213.0 

232.2 



S °4 

51.0 
49.0 
51.0 
98.0 
53.0 
46.0 
46.0 
40.0 
38.0 
51.0 
46.0 
13.0 
46.0 

48. 3 



ca/ 

Mg 

3.03 
5.18 
1.68 
5.05 
3.14 
0.93 
1.30 
3.29 
2.73 
3.07 
1.71 
3.73 
1.95 

2.83 



HCC3/ HC0 3 / 



S0 4 

4.27 
4.03 
4.08 
2.62 
5.25 
3.96 
4.07 
5.93 
5.61 
4.65 
5.70 
16.54 
4.44 



CL 

4.35 
2.39 
4.39 
0.74 
1.54 
1.90 
1.63 
1.39 
8.28 
2.46 
3.02 
6.90 
4.72 



5.47 3.36 



so 4 / 

CL 



1.02 
0.59 
1.08 
0.28 
0.29 
0.48 
0.40 
0.23 
1.48 
0.53 
0.53 
0.42 
1.06 

0.65 



ERROR % 

-1.61 
-0.79 

0.37 

0.02 

0.42 
-1.07 

0.10 
-1.60 
-1.12 
-0.97 

3.62 

1.76 

0.74 



SATC SATD SATG SATM 

0.481 0.490 -1.823 

0.580 0.453 -1.797 

0.489 0.761 -1.897 

0.635 0.581 -1.345 

0.473 0.464 -1.761 

0.621 1.265 -2.379 

1.019 1.873 -2.594 

0.519 0.534 -1.861 

0.416 0.406 -1.993 

0.387 0.298 -1.857 

0.657 1.095 -1.898 

1.229 1.880 -2.412 

0.611 0.939 -1.909 

0.624 0.849 -1.964 -0.254 



CHAMBERS AREA GROUP D SAMPLES 



No. Title 

1 WC8F 

2 DW21 

Mean 



Ca 



Mo 



300.0 75.0 
104.0 77.0 

202.0 76.0 



Na 



305.0 
540.0 



422. 5 



CL 



38.0 1160.0 
68.0 463.0 



53.0 



!11.5 



ALK 


so 4 


Ca/ 

Mg 


HCO3/ 
SO4 


HCO3/ 
CL 


SO4/ 
CL 


ERROR % 


SATC 


SATD 


SATG 


SATM 


121.0 
1140.0 


88.0 
60.0 


2.43 
0.82 


1. 32 
18.24 


0.07 
1.74 


0.06 
0.10 


-2.21 
-1.36 


1.997 
0.847 


3.567 
1.828 


-1.414 
-2.023 


1.091 

0.503 



630.5 74.0 1.62 



9.78 



0.91 



0.08 



1.422 2.698 -1.718 



0.797 



CHAMBERS AREA GROUP C SAMPLES 



No. Title 



WC8R 
DW19 
DW20 
DW24 
DW2 5 
DW26 
DW27 
DW28 



Mean 



Ca 

320.0 
80.0 
65.0 

166.0 
80.0 

109.0 

142.0 
143.0 



Mg 

59.0 
35.0 
25.0 
46.0 
14.0 
16.0 
10.0 
10.0 



138.1 26.1 



Na 

63.0 
18.0 
14.0 
28.0 
13.0 
17.0 
37.0 
43.0 

29.1 



6.7 

5. 3 

1.9 

4.2 

2.0 
1.7 

5.2 

5.4 

4.0 



CL 

132.0 
20.0 
7.0 
53.0 
25.0 
30.0 

63.0 
63.0 

49.1 



ALK 

1040.0 
306.0 
239.0 
326.0 
225.0 
269.0 

385.0 
386.0 

397.0 



so 4 

22.0 

49.0 
50.0 

* 

13.0 

45.0 

28.0 
30.0 

33.8 



CA/ 

Mg 

3.29 
1.39 
1.58 
2.19 
3.47 
4.13 
8.61 
8.67 

4.17 



HCO3/ HCO3/ 



so 4 

45.37 
5.99 
4.59 

* 

5.74 
5.74 

13.20 

12.35 



CL 

5.58 
10.84 
24.19 
4.36 
6.35 
6.35 

4.33 

4.34 



14.84 8.30 



so 4 / 

CL 

0.12 
1.81 
5.27 

* 

1.11 
1.11 

0.33 
0.35 



ERROR % 

-2.54 
-0.79 
-0.53 
25.16 

2.57 

2.59 
-2.04 
-0.74 



SATC 

2.967 
0.464 
0.296 
1.156 
1.315 
0.745 
0.595 
0.597 



SATD 

5.002 
0.800 
0.403 
1.987 
2.081 
0.f'85 
0.275 
0.276 



SATG 

-2.202 
-1.913 
-1.937 

* 

-2.418 
-1.800 
-1.937 

-1.907 



1-34 1.017 1.464 -2.016 

* Sample depleted - analyses not run 



SATM 

1.557 
-0.142 
-0.371 
0.352 
0.287 
-0.339 
-0.798 
-0.799 

-0.032 



WILTON CREEK- MAP 2 
SURFICIAL GEOLOGY 



76°56' 5.5 



76°37' 




/r,"!)f, 



LEGEND 



CENOZOIC 

PLEISTOCENE 
RECENT 



Swamp and alluvial deposits: tho alluvial deposits are limilorl in 
areal extent, being confined to the vicinity of Wilton Creek. Mainly clay, 
silt and fine sand. The swamp deposits, comprised of peat, muck, silt 
and fine sand, are more extensive. 



WISCONSINAN 



Fine-grained glaciolacustrine deposits: mainly clay and silt. These 
doposits are of limited thickness and discontinuous in areal extent. 



Coarse-grained glaciolacustrine deposits: well stralified gravel 
rind sand beach and bar deposits, representing a brief water stand. 



Glacial deposits: primarily sandy-silty till in ground moraine and 
drumlin deposits, 

Glaciofluvial deposits: stratified sand, silt and gravel with extreme 
range in grain sire, sorting and structure. 



SYMBOLS 

Geological boundary, approximate 



300 — Topographic contour 



-^N76"E 



Sand and gravel pit 



Glacial striae 



^$r Drumlins 



_£*=T Drumllnoids 



8 Overburden sample locations 



SOURCES OF INFORMATION 

Surficial geology compiled by G. Funk, 1 973. 

Information supplemented with data from: 

Mirynech, E„ 1 967-68, Unpublished field notes for Bath-Sydenham map area. 
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WILTON CREEK — MAI* 1 
BEDROCK GEOLOGY AND TOPOGRAPHY 



LEGEND 



PALEOZOIC 
ORDOVICIAN 



Verulam Formation 



Bobcaygeon Formation: Upper Momber 



Bobcaygeon Formation: Lower Member 



Gull River Formation: Members B.C.D 



Gull River Formation: Member A, upper submember 



Gull River Formation: Member A, lower and middle submembers 



Shadow Lake Formation 



PRECAMBRIAN 



Precambrian rocks, undifferentiated 



SYMBOLS 



Geological boundary, approximate 



05° Bedrock surface contour, interval 50 feet 



Outcrop exposure 



Quarry 



SOURCES OF INFORMATION 



Bedrock geology and topography as modified from Liberty, 1971, by G. Funk, 1973, 
from cited reference, water-well records assembled by the Ministry of the Environment, 
field surveys and personal communication with B.A. Liberty. 

Liberty, B.A., 1 971 , Paleozoic geology of Wolfe Island, Bath, Sydenham and Gananoque 
map-areas, Ontario; Geol. Surv. Can., Paper 70-35. 
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SAND AND GRAVEL DEPOSITS 



76°37' 







AMHERST -I-^LAND^ 



76°S6' 55 



LEGEND 



Exposed deposits 



Buried deposits 



SYMBOLS 



Geological boundary, approximate 



Geological boundary, assumed 



Sand and gravel pit 



SOURCES OF INFORMATION 



Sand and gravel deposits and their areal extent, as delineated by G. Funk, 1973, from 
water-well records assembled by the Ministry of the Environment and from field surveys. 
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Thickness of Overburden 

0-10 feet 



10-20 feel 



20-30 feet 



30-40 feet 



More Ihan40feet 



SYMBOLS 



Line of equal overburden thickness (isopach), interval 1 feel 



12 Observed overburden thickness 



.20 Overburden thickness from water- well records, in feel 



SOURCES OF INFORMATION 



Overburden thickness compiled by G. Funk, 1973. from field surveys and water-wel 
records on file with the Ministry of the Environment. 
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WILTON CREEK— MAP 6 

INSTRUMENTATION AND 
WATER SAMPLE COLLECTION SITES 



NN 



AR 



GR 



f* 



SG 



RA 



RS 



HT 



SN 



ES 



AT 



DWC-3 



LEGEND 

Hydrometric station, natural control, non-recording 

Hydromotric station, artificial control, recording 

Ground-water well, recording 

Piezometer (asterisk replaced by number of piezometers in nest) 

Snow gauge 

Rain gauge, recording 

Rain gauge, standard 

Hygrothermograph 

Radiometer 

Evaporation station 

Air temperature 

SYMBOLS 

Water sample collected from MOE observation well 
Water sample collected from domestic well 



ST-12 Water sample collected from stream 



SOURCES OF INFORMATION 

WBter samples collected on field surveys in 1969 and 1 973. 
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Figure 3. Gamma logs for Ministry of the Environment observation wells on a cross-section near Violet. 
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